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This serves as a final report for grant NAG5-1554 entitled Energization and Transport of 
Ions of Ionospheric Origin in the Terrestrial Magnetosphere. The work has been predominantly 
focused on ion outflows identified in two data sets: 1) Prognoz 7 and 2) Dynamics Explorer. 

The Prognoz analysis has been carried out in conjunction with Dr. Oleg Vaisberg of IKI. 
The study analyzed ion densities, temperatures, and flow velocities in the magnetolail. The work 
performed under this contract consisted of developing a program to load the raw data, compute 
the background subtraction of a strong sun pulse, and use the net counts to calculate the low 
order moments of the distribution function. The study confirms the results of ISEE that the the 
cusp is a major source of plasmashect plasma and goes beyond this to discuss the use of ion 
velocities as a way to examine the motions of the magnetotail. The initial work will be reported 
in Advances in Space research and is included in abstract form as Appendix A. In addition. 
Appendix A also contains a draft of a paper by Delcourt et al. which will be submitted to 
Geophysical Research Letters in June. 

The work on the DE/Relarding Ion Mass Spectrometer data set represents the bulk of the 
effort. The work carried out under this grant is separated into two categories: 1) classification 
of low-energy ion flows from high-latitudes, and 2) studies of the polar wind. The work on 
classification of ion flows has consisted mainly of support work on studies that begin before my 
departure from NASA Marshall Space Flight Center. The major publications that have resulted 
from this work are: 1) the determination of plasma fluxes from the cleft ion fountain Pollack el 
al. [1990] and the recent thorough analysis and characterization of DE/RIMS high-latitude ion 
outflows by Giles et al. [1995]. These reprints are included as Appendix B of this report. 

The polar wind is indeed a component of the low energy ion outflow at high latitude. 
However, the extensive theoretical work on this subject puts it in a category all to itself. The 
analysis of the polar wind by the RIMS group in the late I980’s culminated in the work of 
Chandler et al. [1991] (see Appendix C). This work reports the flux variability as a function of 
season, magnetic activity, etc. and was partially supported by this grant. In addition, much effort 
has been expended under this grant to complete a follow on study of the thermal structure of the 
polar wind. Extensive display tools and analysis software have been developed and used in an 
attempt to carry out this thermal analysis. However, the difficulties of performing multi- 
parameter, non-linear least squares fits has so far proved to be insurmountable. Therefore, the 
present work (see notes in Appendix C on work still in progress) uses a constrained fit scheme 
that combines the ion densities and flow velocities derived from Chandler et al. [1991] and a 
spacecraft potential derived from an empirical relation to the total ion density (utilizing the 
complete RIMS data set) to determine the remaining fit parameter, the ion temperature, via a 
least squares fit to the RIMS data. Reports of this work are anticipated in the near future. 
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ABSTRACT 

The plasma spectrometer SKS on the Prognoz 8 satellite was able to measure the very dilute 
plasma in the tail lobe at distances 15-20 Earth radii about 50% of the time. The energy of the 
ions, being in the range of several hundred electron volts, decreases as the satellite moves from 
the high-latitude magnetopause to the plasma sheet. The velocity spread of the ions is narrow 
in the direction of motion, corresponding to a temperature of several eV. The temperature in the 
transverse direction is 2 to 3 times higher. Comparison with magnetospheric convection models 
suggests that these ions originate at low altitudes in the region near the cusp and arc separated 
in the magnetospheric velocity spectrometer. Observed velocity variations along the spacecraft 
trajectory suggest the existence of a characteristic velocity profile and a characteristic 
temporal/spatial scale. Comparison of simultaneous E/q measurements with mass-selective 
detector measurements of protons and the alpha particles suggests that lobe plasma at these 
distances consists primarily of oxygen ions. The average ion flux density is about 3X10 5 cm' 2 
s' 1 , giving estimated total flux over the lobe about 10 25 s' 1 , comparable to cusp/cleft source. 
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SIGNATURES OF IMPULSIVE CONVECTION 
IN THE MAGNEi OSPHERIC LOBES 


D. C. DeJcourt J , J A Sauvaud* , (). L. Vaisberg* ,I,A. Avanov\ 
.1 I. Burch* , and J H Waite, Jr 1 2 * 4 


Abstract. Wc examine the injections of tailward flowing 0 + ions measured by the 
l‘llOGNOZ-8 spacecraft in the magnetospheric lobes using three-dimensional single-pai tide 
codes. The repeated energy-lime injection patterns observed are unexpectedly well correlated 
with AE variation* during a time interval of about 7 hours. Such a correlation leads us to 
consider the effects of transient intensification of the large-scale magnetospheric plasma 
convection. Using a simple model of time-varying electric field, we show that the source of the 
ions lies near ihe day side cusp. This source extends over a narrow (a few degrees) latitudinal 
interval in the vicinity of the last closed field line, likely involving ions expelled into the 
magnetosphere after being transversely accelerated We demonstrate that, downstream of the 
ejection sito, impulsive convection leads to prominent structuring of the tailward streaming 
particles with abrupt energy variations on small length scales. We sliess that the changes in 
large-scale plasma convection are linked (u processes internal to the magnetosphere occurring 
during substorms. The agieeinent between the numerical results and the data provides further 
evidences that there exists a close relationship between the AE magnitude and the convention 
electric field. 
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1. Introduction 

On August 27 1981, the PROGNOZ-8 spacecraft traveling inbound in the magnctospheric 
lobe (northern hemisphere) recorded several injections of ions with energies in the 200-eV to 
1-keV range and propagating in the antisolar direction. In a previous study, Vaisberg ei cd. 
[tyyb] identified these tailward flowing ions as cold ion beams mainly consisting of 
ionospheric (V and suggested they originate from the cusp region. The bottom panel ol'l’lale 1 
shows Ihe energy-time spectrogram recorded during this pass. In this panel, it can he seen that 
the PROGNOZ-8 spacecraft enters into the magnetosphere at approximately 06:30 UT as 
identified by the intense flux of energetic tailward flowing particles in the plasma mantle 
adjoining the magnetopause [Rosenbaucr ei al., 1975], Subsequently, as the spacecraft travels 
through the lobes in the vicinity of the midnight meridian plane, repeated structures arc 
noticeable which are well separated in time and nearly shaped as inverted Vs. These structures 
occur from -12.00 UT (Z ~ 18 Rjt) until '-23,00 UT (Z - 8 Rj?) Thu spacecraft then 
encounters the plasma sheet and ultimately the radiation belts. 

It is of interest to note that the occurrence pattern of the ion injection structures and the 
simultaneous Alt, variations shown in the top panel of Plate I are quite similar in shape As a 
matter of fact, a cross -correlation analysis between the AK magnitude and the energy of the 
maximum 0 + flux in the spectra reveals a correlation coefficient of -0.7 between 13:00 and 
19:00 UT with a timo lag smaller than - 10 minutes (the AK change proceeding the energy 
change). This correlation suggests that tire injection structures arc linked to global changes in 
the magnetosphere and intimately related to subslorm processes. As these AK variations aic 
indicators of substoruis which lead to tenipoiuiy intensification of The large-scale 
magnetospheric plasma convcciiuu, the purpose of this paper is to revisit more quantitatively 
the results of Vaisberg et al. (19951 usinft three-dimensional calculations of single-particle 
trajectories in presence of a time-varying electric field 

2. Origin of Tailward Plowing t) + 

In view of the good agreement between the A F. variations and the observed injection 
pattern, we start with the assumption that the cross-tail potential drop can be directly related to 


1 
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the AE magnitude. Such an assumption dearly is oversimplified. In particular, large changes in 
the AE magnitude provide evidences of modifications in the magnetospheric current system, 
and the above assumption does net cxplicitcly take into account the induced electric field 
which should result from concurrent dynamical reconfiguration of the magnetosphere. Still, we 
will see from the simulation results that, in first approximation, tire above assumption is not 
unreasonable in the magnetospheric lobes and that it leads to injection signatures which are 
consistent with those observed, As Illustrated in the two upper panels of Figure 1, we focused 
on the 1U:UU to iy UO UT sequence of the PROGNOZ-8 pass and modeled the cross-tail 
polenlial drop (<J») in a nearly linear manner with the AE magnitude, imposing threshold values 
of 40 kV for Alt £ 200 nT and 200 kV for AE i 1000 nT. The variations thus obtained arc in 
qualitative agreement with previous detailed analysis of the relationship between the 
magnetosphoric electric field and the current in the ionosphere [e g., Kamido and Baumjohann, 
1985; Richmond ct al, 1990; G. Lu, private communication, 1995], Also, the resulting three- 
hour averaged potential drops (heavy dashed line in the lop panel of Figure 1) are comparable 
to those derived using th t Maynard and Chen [1975] equation and appropriate values uf the 
Kp parameter, viz., Kp = 3 A from 9 00 to 12:00 UT, Kp - 5 + from 12.00 to 18 00 UT. and Kp 
= 4" 1 from 1 8:00 to 21 :00 UT [Coffey, 1981], 

In order to perform single-particle analysis of the energy-time spectrogram in Plate i , we 
sorted out the energy of the maximum G + flux during each 4-minute scan of one of the 
PROGNOZ-8 spectrometers. Note that the three spectrometers onboard PROGNOZ-8 were 
looking within ~ +/-7.5 0 from the Sun direction (i.e., close to tho magnetic field direction). The 
energy profile is very similar for all three instruments and is illustrated in the third panel from 
top in Figure 1 (heavy solid line). The repeated invcrtcd-V like structures mentioned above arc 
clearly apparent in this panel, and their modulation according to the AE magnitude is readily 
noticeable by comparison will) the second panel fiuui lop. To identify the origin of these 
lail ward flowing ions, we performed numerical trajectory calculations backward in time from 
the spacecraft location In the Isyganenko [1987] model (referred to as T-87 hereinafter) fcach 
sample of the observed distribution was assigned a test O’* with corresponding energy, 
assuming the particle to be nearly aligned with the magnetic held (i.e , nearly zero magnetic 
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moment) in the guiding center reference frame. Under such conditions, the particle pitch angle 
in the Earth's reference frame is essentially due to the local instantaneous ExB drift velocity. 
These pitch angles were found to be comprised between 5° and 15° at the spacecraft altitude 
(fourth panel from top in Figure 1), that is, in a range consistent with the direction of view of 
the FROGNOZ-8 detectors ( Vaisberg d al. [1995]). As for the location, it must be stressed 
that the spacecraft encounter with the magnetopause at ~6;30 UT occurs somewhat lower in Z 
height than the position obtained in the T-87 model (sec open square in the left panel of Figure 
2). This difference is of importance when tracking the origin of the tailward flowing panicles 
as it can load to unrealistic mapping into the ionosphere. In the numerical calculations, we 
adjusted the model magnetopause and Ihe observed one by heightening the spacecraft position 
by about 1.5 in Z height (sue left panel of Figure 2), that is in practice, by retarding the 
spacecraft orbit by about 3.8 hours. 

The results of these calculations are shown in Plate 2 which shows the computed O* 
trajectories (color-coded according to their arrival lime al the spacecraft.) until the first mirror 
point encountered during backward tracing. These mirror points were all found below 1 if/,; 
altitude. Also, the tluue lowci panels uf Figure 1 present the magnetic local lime (MET), die 
invariant latitude (1LAT) and the time of flight of the particles at mirror point. It can be seen in 
these panels that most of the ions originate from a narrow region in the ionosphere, both in 
Ml.T (between -10 00 and -I I 00 hours) and in TLAT (between -70° and -80°) Here, it 
should be pointed out that the numerical results are fairly sensitive to the particle pitch angle at 
the satellite position as well as to the model adopted for the time-varying electric field. Trial 
and orror computations showed that only particles nearly aligned with the magnetic field and 
drifting in an electric field nearly proportional to the AE magnitude originate from the topside 
ionosphere during the time interval under consideration. Other combinations gave unrealistic 
sources of O 4 ' located well above the ionosphere or even in (.lie oulci magnetosphere. 

On the othei hand, the vertical bars in the third panel from top reveal initial energies for 
these O 4 between -100 eV and several hundreds of cV These results clearly corroborate the 
suggestion made by Vaisberg ct al. [19951 that these ions originate from the cusp region. In 
particular, the T-87 model used In these calculations {Kp = 4.4 + ) places the last field line which 

4 
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is dosed in the dayside sector near 75° ILAT, and it can be seen that the ions detected by 
PROGNOZ-8 originate fiom within a few degrees on eilhej side of this last dosed field line. 
Significant wave activity is known to occur in this region of the high-latitude ionosphere fe.g., 

^ J,undin, 1988] and from the present results (i.e.. characteristic energy of the order of a few 

j 

j hundreds of ev and pitch angle of 90° below 1 Kj,; altitude), ion conics locally produced via 

transverse healing are iikely to be the source of the observed tatlward flowing ( ) + 

i 

L 

! 3. Impulsive Magnetosphere Convection 

] An interesting feature of the PROGNOZ-8 data is the sharp inverted V-like structuring of 

[ the ion energy (third panel from top in Figure 1). In particular, clear interruptions arc 

noticeable between two successive structures which cannot be interpreted simply from 
; geomagnetic mass spectrometer [)x>ckwood et al , 1985] effects in steady stale. Indeed, ExB 

drift in presence of a constant dawiwo-dusk potential diop leads to spatial dispersion of 

i 

pai tides with diflinunl uuiagies (or masses) but in a continuous manner. This Clearly is 
incompatible with the dispersion pattern observed. To account for the abrupt variations 
portrayed in Figure 1, one must consider temporal effects within the geomagnetic mass 
i spectrometer I o do so, we performed numerical trajectory calculations forward in lime from a 

\ given source in the topside ionosphere, adopting a time-varying potential drop as shown in the 

1 top panel of Figure 1 . In view of the two bottom panels of Figure 1 , the source in the 

ionosphere was taken to extend over ILAT on either side of the last closed field line at 
j 1 1 .00 MLT. Test 0 + were launched from this line source at distinct times (by steps uf 10- 

• minutes) throughout the 10:00-19:00 UT interval. 

Examples of the trajectories achieved aic given in Figure 2. In this figure, one of the 0 + 

i 

• ( s0 'id line) is initialized at 76° ILAT while the other (dashed line) is launched 10 minutes later 

1 fiom » higher invariant latitude (77° ILAT). As a general comment, it can be seen in the upper 

ri Hht panel of Figure 2 that, in both cases, a significant energy gam (from 300 eV up to -600 
j tV) 1$ obtained filler about 10-minutc time oi flight bs the panicles travel through a region of 

enhanced Held line curvature, As can he seen by comparison with the center right panel, this 

) 

energy gain essentially occurs in the parallel direction as a result of curvature drift through the 
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convection electric field equipotentials [e g., Cladis , 1986]. It is also apparent that phasing 
between the panicle orbits and the time-varying convection eleouic field (lower light panel of 
Figure 2) leads to substantial variations in the net par ticle path Most notably, the 0 + launched 
from a higher 1LAT (dashed line) which should intercept the spacecraft at a later time and a 
lower Z height in the magnctosphcric lobe, is actually recorded at a position similar to that of 
the O* originating from a lower II A T (solid line). This phasing effect combined with the ahove 
curvature-related acceleration can provide the mechanism for the well-separated inverted V- 
likc structures in Figure 1 . 

This can bo appreciated in Figure 3 which shows the energy variations obtained along the 
satellite pass assuming either a constant (100 kV) dawn to dusk potential drop (top panel) or a 
time-varying one (second panel from top). For simplicity, a 300-cV initial energy was 
considered in both cases. In the steady-state case (top panel), a continuous pattern is obtained, 
exhibiting a monotonous energy decrease from ~700 eV ai 10.30 UT duwn u> ~450 eV at 
13:30 UT. 'these variations result ftom enhanced Add line curvature in the outer cusp region, 
and thus larger energization for particles launched from lower lLATs As mentioned above, 
this pattern is incompatible with that observed (light lines) even though the average energy is 
fairly well reproduced In contrast, in the second panel from top, narrow 0 + structures with 
increasing and decreasing energies are noticeable from -13:00 to -18 00 UT, which resemble 
those observed. In particular, the energy modulation follows quite well the AK variation. This 
agreement is essentially qualitative as the amplitude of the energy change is substantially 
smaller than that observed. This latter effect may be due to the mono-cncrgetic character of the 
model source considered. As for the time interval before -13:00 UT, the computed pattern 
sharply differs from the observations. However, we note that litis period of lime corresponds 
to a fairly weaker magnetic activity (Kp - 3 + ) which may lead to distinct characteristics of the 
uulfluw [c.g., Yau e( ul . , 1986), 

Returning to the afternoon sequence, it must be stressed that the good agreement was 
obtained by shifting the numerical results by -13 minutes toward earlier times or, equivalently, 
by considering that the injections respond to variations in the convection rate (Ah' magnitude) 
with a time delay of the order of 13 minutes. Within the limits ofihe model, this estimate is 
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compatible with that obtained by direct comparison of the energy spectra and the A/i variation 
profile (sec introduction). This time delay is significantly smaller than the time of flight of the 
particles from the ionosphere to the spacecraft (of the order of 40 minutes; see Figure 2), but it 
is comparable to the time required for an 0 + to travel over the distance between the spacecraft 
and the legion of enhanced field line curvature where energy can be gained. In other words, 
what the numerical model suggests is that (he stiuctuied injections observed do not result front 
a temporal modulation of the O 1 source in the cusp region (which would take too long to 
propagate down to PROGNOZ-8) but from the time-varying electric field acting upon the 
particles during their magnetospheric transport. This mechanism which operates in the 
magnetophone lobes is reminiscent of that discussed hy Mauk 11986] in the near- Hart h plasma 
sheet whereby' phasing between a surging electric field and the particle longitudinal motion 
leads to repeated injection patterns. 

Finally it should be pointed out that, due to the strong variability of the particle drift paths 
with time (as illustrated at left in Figure 2), temporary flux dropouts arc obtained in Figure 3 
which arc consistent with those featured in the data. These flux dropouts coll for a limited 
extension of the source region in the ionosphere. Indeed, it can be seen in the second panel 
Hum bottom in Figure 3 Llutl each injection involves ions originating from lower ILATs owing 
to the intensifying magnetospheric convection, in contrast, after each injection, the upper 
latitudinal limit set on the source region (which extends from 71° TLAT up to /y v ii.aT) 
prohibits particle transport from higher II. A I s to the spacecraft. Note moreover in the bottom 
panel of Figure A that the model t)* have pitch angles in the range of the PROGNOZ-8 
instrument (i.e., within -15° from the magnetic field direction). 

4. Summary 

The intent of this paper was to provide quantitative insights into the tailward injections of 
O 4- ions reported by Vaisberg et al. [1995] using numerical trajectory calculations. The clew 
correlation we found between the observed injections and die magnetic activity (as viewed 
front the AE variation) suggests to consider the effect of sudden enhancements in the 
magneiosplieiic convection rate. The interpretation framework developed in this study clearly 
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'is a swnpftftvd description of panicle transport during such events. Nevertheless, there aic dear 
indications that impulsive convection is at work and is responsible for the observed injections, 
namely: (1) numerical simulations of this effect place the source of the tailward propagating 
ions in the vicinily of the daysidc cusp, i.c „ in a region known to be an active source of Cf. 
and (2) steady-state plasma convection cannot Jead to a dispersion pattern compatible with the 
observations. In contrast, the repeated injection structures as well as the flux dropouts between 
them are well reproduced if one considers large electric Held variations on short lime scales. 
These variations lead to distinct energy gains during transport from tho ionosphere as well as 
inlricalcd drift paths which occasionally intercept the spacecraft It may be that part of this 
electric field intensification is of inductive nature though the present study lumps it with the 
large-scale convection. These results cleaily exemplify the role of subsluitu expansion 1 dated 
electric fields in the iuagneiuspheiic lobes 
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Plato Captious 

riato 1. (Top) AE index versus time, (Bottom) Color-coded energy-time spectrogram 
recorded by the PROGNOZ-8 spacecraft on August 27 19H1 This spectrogram relate to 
panicles propagating tailward with pitch angles smaller than 15°. 


Plate 2. Model O" trajectories in the noon-midnight meridian plane. These trajectories are 
color-coded according to the arrival time at the PROGNOZ 8 location (heavy solid line) 
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A Survey of Upwelling Ion Event Characteristics 

C. J. Pollock . 1 M. O. Chandler . 1 T. E. Moore.' J. H. Waite. Jr . 2 
C. R. Chappell , 1 and D. A. Gurnett' 


Ionospheric ion upwelling in ihe vicinity ot the davsidc cleft has been studied, based primarily on data from 
the Dynamics Explorer I spacecraft. Using retarding ion mass spectrometer low-energy ion data and plasma 
wave instrument dc electric field data, bulk ion plasma parameters, including ion species density and field- 
aligned bulk velocity and tlux. have been derived at points within a number of observed upwelling ion events tor 
the ion species H * . He * . O * . and O * V The ion species bulk parameters near the source latitude are examined 
and compared. Wc find that the upwelling plasma is rich in O * . which typically comprises of the particle 

density, followed by H at somewhat less than I09K and then He* and O* *. each comprising —I# of 
the upwelling ion panicle density. The upwelling O* ion tlux is also commonly dominant over that associated 
with the other species, with normalized values near the source region which are typically near 10 s * cm : s 1 . The 
fractional upward H tlux is not as small as the fractional H* density due to the much larger H* upward 
flow velocities. Integration ot the product of the normalized upward ion species tlux and the upwelling ton 
occurrence probability (Lockwood et al.. 1985) over the source area yields an estimate of the source strength of 
this low-altitude cleft region magnetosphenc plasma source of 2.6 x I0 :fi ions s 


I . Introduction 

During recent years, it has been consistently demonstrated that 
the Earth's ionosphere acts as a significant source of magneto- 
spheric plasma [Shelley et al.. 1972; Young et al 1977; 
Ghielmetri et al.. 1978; Johnson. 1979; Balsiger et al.. 1980; 
Gorney et al.. 1981 ; Lundin et a/.. 1982; Collin et al., 1984; You 
et al.. 1984, 1985; Lennartsson and Shelley, 1986). In fact, it has 
recently been argued f Chappell et al.. 1987) that ionospheric 
plasma may fully populate the magnetosphere, with no significant 
solar wind source being required to account for observed mag- 
netosphenc charged panicle populations. Recognition of the 
importance ot the ionosphere as a source of magnetosphenc 
plasma began with observations of heavy ions, previously 
believed to be gravitationally bound to low geocentnc altitudes, at 
magnetosphenc altitudes. As such observations continued to be 
reported, in diverse magnetosphenc regions such as the ring 
current [Johnson et al.. 1977). the plasma sheet {Peterson et al., 
1981), and the plasma sheet boundary layer and tail lobes 
| Eastman et al., 1984|. it became apparent that escaping plasma 
transport from very low altitudes must be quite common. At the 
beginning of the 1980s, however, the questions of transport 
mechanisms and specific source regions were only beginning to 
be observattonally addressed. 

Observations of outflowing ionospheric ions have been 
reported by a number of authors in the last decade. Outflowing ion 
distributions which have apparently evolved, beginning with 
transverse (to B) acceleration and followed by magnetic mirror 
folding, have been observed near I Rfr [Sharp et al.. 1977) in the 
topside ionosphere [ Klumpar, 1979) and at sounding rocket alti- 
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tudes [Whalen et al.. 1978; Yau et al.. 1983; Moore et al.. 1986a; 
Kintner et al.. 1989). Since the launch of Dynamics Explorer I 
(DE I), observation of low-energy ionospheric plasma outflows 
has become very common. The realization that the daystde mag- 
netospheric cleft region represents a significant source of outflow- 
ing plasma began early in the DE mission, based on backward 
trajectory analysis of upflowing beams observed in high-altitude 
(2-4 R f ) polar cap field lines ( Waite et al.. 1985). This realization 
reached matunty with the work of Lockwood et al. ( 1985) who 
used DE 1 retarding ion mass spectrometer (RIMS) data to study 
the statistical occurrence of various types of ionospheric outflows. 
That work identified a persistent (occurrence probability exceeds 
0.5 in source region) source of plasma emanating from the day- 
side auroral zone, showed that this source apparentlv favors the 
prenoon region of magnetic local time, and showed that there is a 
strong relationship between the low-latitude edge of the source 
region and the location of (usually downward) field-aligned 
currents. They labeled the observation of the distinctive RIMS 
spin-time signature associated with this outflow as an upwelling 
ion event (UWI). Moore et al. [ 1 986/?) and Waite et al. [1986] 
presented case studies of an UWI, observed near 1930 UT on day 
71 of 1982. They derived bulk flow parameters for several species 
and demonstrated that the upwelling ion plasma core ts hot and 
conically shaped, having undergone transverse energization at 
some lower altitude. They also noted a strong shear in the 
transverse electric field associated with the event. 

The purpose of the present study is to more broadly investigate 
the bulk properties of UWls, using RIMS and plasma wave instru- 
ment (PWI) quasi-static electric field data. The results of the study 
yield quantification of the densities and the upward field-aligned 
velocities and fluxes, for the ion species H + , He*, O". and 
O" + , associated with these flows. These parameters, as measured 
near the lower-latitude edge of the events, where the flows arc 
most intense, will be presented. All of the available altitude 
normalized ion tlux measurements from within the upwelling ion 
events are combined with the statistical occurrence probabilities 
derived by Lockwood et al. (1985) to produce spatial source 
strength distributions. Integration over these distributions yields 
an improved measure of the total low-energy ion outflow rate 
from this source region. 
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2. Instrumentation 
2.1. The Spacecraft 

The observations presented in this report were gathered using 
instruments flown aboard the Dynamics Explorer I I DE I ) space- 
craft. DE 1 was launched on August 3, 1981 . into a highly ellipti- 
cal polar orbit with apogee and perigee altitudes at 23.300 km 
*4 06 Ry . geocentric I and 570 km. respectively. The orientation 
ot the orbit plane is fixed in the geocentric equatorial inertial 
iCjEI) frame of reference, so that in the geocentric solar ecliptic 
iGSE) frame of reference, it is seen to rotate westward at a rate of 
0.98 u /dav. This provides for approximately one full 360° rotation 
of the orbit plane in geocentric local time per year. Further, the 
DE l line of orbital apsides rotates about the negative orbit plane 
normal at a rate ot l20°/year. so that the geographic latitudes of 
apogee and perigee drift retrograde at that rate. The spacecraft 
orbit, then, very nearly repeats itself, in the Earth's frame of refer- 
ence. once every 3 years. These aspects of the spacecraft orbit 
influence our ability to study UWIs. not only because the source 
region has a limited geocentric solid angular extent, but also 
because the flowing plasma seems to evolve substantially with 
ulntude. Therefore, the plasma characteristics which give rise to 
die distinctive RIMS UWI signature appear to exist only in a small 
region of space sampled by the spacecraft for several months once 
every 3 years. Furthermore, the orbit forces a phase-locked 
altitude-local time sampling which may bias the occurrence statis- 
tics ol the local time extent ot the upwelling ion source region. 
Other aspects of the DE mission have been described elsewhere 
[Hoffmann and Schmerling , 1981 (. 

2.2 Retarding Ion Mass Spectrometer 

Positive ion parameters presented in this report were obtained 
using the DE 1 retarding ion mass spectrometer < RIMS) [Chappell 
ct a! .. 1981 1. The RIMS instrument consists of three similar 
sensor heads, one each viewing parallel and antiparailcl to the 
spacecraft spin axis (negative and positive Z heads, respectively) 
and the third viewing in the spin plane (radial head). The two Z 
heads feature circular fields of view with half angles of 45°. while 
the radial head views a field which has a 45'"' half angle in the 
plane perpendicular to the spin plane and a 20° half angle in the 
s P m piane. Since the DE 1 spin axis is nearly perpendicular to 
the geomagnetic field, the fields of view of the two Z heads arc 
centered nearly perpendicular to that field, and the field of view of 
the radial head sweeps through a nearly full range of pitch angles 
with the spacecraft spin. All three sensor heads were designed to 
discriminate among incoming positive ions on the basis of energy 
per charge (0-50 cMu/), using a retarding potential anaivzer 
( RPA). and on the basis of mass per charge 1 1-32 aniu/i/). using a 
magnetic mass spectrometer. Further, each sensor utilizes two 
detectors that are sensitive to ions with mass ratios of 4 to 1 such 
dial, for example, while one detector is viewing O * . the other is 
\ icw ing He * . The RPA on the radial head failed in late 1981 . with 
the retarding grid having become electrically grounded to the 
spacecraft chassis so that, subsequent to the failure date, all ions 
with energies per charge greater than any positive spacecraft 
potential w'ere passed into the magnetic mass spectrometer. The 
mass spectrometer, in turn, is characterized by a finite energy 
passband. whose width is dependent upon the species mass per 
charge and vanes from approximately 250 eV for H * down to 16 
vY for O " The upwelling ion densities, field-aliened bulk 
velocities, and fluxes presented below were derived using ion data 
trom the RIMS radial head. 


2.3. Quast-Siiitn Electric f ields 

Quasi-static electric field data presented in this report and used 
in deriving ton nulk parameters were obtained using the / axis 
antenna of the plasma wave and quasi-static electric held instru- 
ments [Shawnan et id., 19X1 |. liiis tubular, (tip to tip) 
antenna is sensitive to the electric tield component perpendicular 
to the spacocrait spin plane. I he geomagnetic Held is nearlv 
(within 5 m l *5 l r ol our eases and alwavs within 12 for our 
data set) contained within the DE I spin plane, so that the / 
component electric field is associated with convective plasma 
drifts which are approximately within the spin piane and perpen- 
dicular to the ambient geomagnetic field. This electric field 
component is prov ided at a rate of 16 samples s 1 over a nominal 
range of 0.5 m\ m 1 lo 2 V m 1 Characteristics of this instrument 
which are pertinent to the analysis techniques employed in this 
study w'ill be described below. 

3. Data Analysis 

3.1. Data Selection 

The observations presented here are taken from a subset of up- 
welling ion events identified by M. Lockwood through a survey 
of the RIMS data as displayed in the DE I microfiche (with format 
similar to Figure 2 a displaying data for the species H ‘ and He " ) in 
conjunction with his study of the occurrence of various types of 
terrestrial ion outflows [Lockwood et al .. 1985). Lockwood et al. 
used a qualitative visual signature identification technique to 
identify these upwelling ion events. The RIMS spin-time up- 
welling ion event signature has been described elsewhere (see. for 
example. Plate 3 of Lockwood et al. |I985|). This signature is 
again illustrated in the O ' data presented in Figure 2a. where the 
count rate in the O ’ channel is shown plotted in grey scale versus 
UT (along with certain orbital parameters) along the abscissa and 
the RIMS radial head viewing direction, w nh respect to the space- 
craft ram direction, along the ordinate. The solid and dashed lines 
running through the spectrogram from left to right indicate view- 
ing angles along the positive and negative geomagnetic field 
directions, respectively. For northern hemisphere passes, ions 
detected while viewing along the positive magnetic field direction 
are moving in the anntield direction and. theretore. awav from the 
Earth. The signature of the event is evident in Figure 2. beginning 
near 2320 UT and continuing until just past 2327 UT. The ramp- 
like upward flux intensification, followed by a sharp cutoff and 
return to ram-dominated flow into the instrument at the low- 
latitude edge of the event, typifies the signature of these events 
used by Lockwood et al. in their fiehc survey. 

Having surveyed the available data obtained between October 
1981 and October 1983, Lockwood et al. identified 86 upwelling 
ion events. These events were identifiable in the tiche for periods 
ot time ranging trom I to 15 min. providing a total of nearly 450 
min of observation within the events. Of these 86 events, we have 
selected 39 for inclusion in the analysis for this report. This selec- 
tion was based, in part, on the requirement that the RIMS instru- 
ment be operating in a mode without aperture bias applied to the 
instrument entrance and on the further requirement that the meas- 
ured Z component of the electric field was deemed reliable. The 
RIMS instrument features an aperture bias capability, w hereby a 
negative electrostatic potential may be applied to a grid at the 
entrance aperture, in order to overcome the effects of a positive 
spacecraft potential, allow ing access of very low -energy ions to 
the instrument. Of the 86 upwelling ion events identified in the 
Lockwood et ai study. 20 were from times when RIMS was 
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operating in the aperture bias mode, and analysis of these events 
has been deferred pending the development of routine analvsis 
loois for this mode. For another 8 events, either the digital RIMS 
data or the attitude and orbit data were unavailable. Of the remain- 
ing 58 events. 16 were rejected because the Z antenna dc electric 
held data were deemed unreliable. This antenna has a length of 
9 m tip to tip. as compared with the spin plane antenna which is 
200 m long. It is known to be subject to stray fields that originate 
trom dittercntial spacecraft charging and can have large amplitude 
relative to the fields ot interest. This left 42 events, from which 2 
were rejected due to a combination of no electric field data avail- 
able and highly variable, difficult to interpret ion data and I was 
rejected due to the apparent presence ot a telemetry dropout dur- 
ing the event. Thus. 39 of the original 86 events, or 45 c /c. have 
been processed to obtain ion bulk parameters along the spacecraft 
track. Of these. 5 have no electric field data, and values of £ ; = 0 
have been assumed. These 39 events, consisting of a total of 
nearly 200 min of observation, constitute the data set used in the 
present analysis. 


5.2. Derivation oj ion Bulk Parameters 

Due to the tailure ot the RIMS radial head RPA. ion enerev 
spectra are not available in the spin plane data. This situation has 
led to the development of an alternative technique for deriving ion 
hulk parameters from the RIMS radial heat data. This geometrical 
technique, developed by Chandler and Chappell 1 1 986 ] . provides 
a measure of the field-aligned ion bulk velocity and field-aligned 
ion flux of all ions of a given species which fall within the energy 
passband of the magnetic mass spectrometer. Derivation of the 
field-aligned bulk velocity depends upon the presumption that 
there are three distinct contributors to the species bulk velocity in 
the spacecraft spin plane and frame of reference (see Figure I ). 
These contributors are the spacecraft-induced ram velocity (V r ). 
F x B drift (V Ex b) and field-aligned drift (V||) velocities. The 
spacecraft ram component is obtained from the orbit and attitude 
data base, while the E x B drift component is derived, using meas- 
ured electric and magnetic fields, as 


V 


K\B 


cE x B 

|B| 2 


I O ' 2 ^ (y) 

|B| s 


(I) 


where, on the extreme right-hand side. E. is the measured Z 
component of the electric field, expressed in mV m and B is 
expressed in gauss. Here, we have used the approximation that the 
magnetic field lies in the spin plane of the DE I spacecraft. As 
was mentioned above, this approximation is always good to 
within 12° for the data used in this study. Then, the total spin 
plane ion species bulk velocity in the spacecraft frame of refer- 
ence is given as 

A / 

vpm pi. mo ““ ( ' rx * id X + ( l' r ,. + y (2) 


where r it is the field-aligned bulk velocity wc arc trying to deter- 
mine. r rr and vvy are the X and Y components of the ion ram 
\elocit\. and 

x = B y - z x x 7 . = -uj s , (3) 

Note that positive values ot ticld-aligned bulk velocity correspond 
to upward and downward drifts in the southern and northern 
hemispheres, respectively. All field-aligned drifts presented in 
this paper are upward, however, and wc have plotted these as 
positive. Here. cu^. is parallel to the DE 1 spin angular 
momentum. The unknown parallel drift velocity is given as 
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Fig. I. The coordinate system and geometry used for computing ion 
species upward bulk velocity from RIMS spin curves and PWI clecinc field 
data is displayed in the top panel. At the bottom, a typical RIMS count rate 
spin curve is shown, and the meaning ot angles identified in the sketch above 
is illustrated. The magnetic pitch angle at which counts arc accumulated is 
displayed as the triangular wave in the bottom panel. 


The angle cx is measured as the difference 

a = ct) p -<i)B (5) 

in spin phase angles between the location (6 P ) of the centroid of 
the distribution of ion counts with spin phase and the direction 
(d>B) along which the positive magnetic field direction is viewed. 
Spin phase angle (6) is that measured about + z. from the space- 
craft velocity vector to the RIMS radial head look direction. In 
order to derive the field-aligned ion flux for a given species, the 
product ot the instrument count rate and calibrated instrument 
sensitivity factors arc cosine weighted, with respect to the spin 
plane velocity vector (given in equation (2). above), and inte- 
grated over spin phase angles, yielding the tlux (J T ) along the 
direction ot the total spin plane velocity vector. The field-aligned 
flux is then obtained as the product 


h = h ■ L— 

Ospin plane! 

and the species panicle density is estimated as 


J T 


_ A 

v ’ll 


( 6 ) 


(7) 


I'Npm plane! 

where iv sp , n p , un< .j is determined using equation (2). having deter- 
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mined v KxB using measured fields as in equation ( I ). the ram 
components using orbit and attitude data in conjunction with mag- 
netic field data, and vji using equation (4) as described above. Our 
estimate tor the species particle density is actually a lower bound, 
which disregards the possible presence of ions at energies lower 
than any positive spacecraft potential. This technique has been 
applied to data obtained in and near the 39 upwelling ion event 
intervals included in this study. We have used I -min data aver- 
ages throughout in analyzing the PWI electric field and the RIMS 
ion data. 

v 3 Error Analysis 

l he above analysis involves tw f o potential sources of error. 
First, lor angles (a) near i t or 0. corresponding to highly field- 
aligned or anti-field-aligned Hows, the analysis produces values 
ot \*u which are unrealistically large and positive or negative, 
respectively, for any nonzero difference (v 0 , - |i‘e*b|). Wc have 
discarded results which yield a total velocity so large that the 
corresponding ion kinetic energy per charge would be too large to 
Ik- passed by the magnetic mass spectrometer. The accuracy with 
a Inch we can determine a is on the order of ±2°. so that 
• tie analysis ot flows which arc in actuality within —2° of being 
lieid-ahgned would be expected to produce values of \' tl whose 
magnitudes are arbitrarily large. We have propagated an uncer- 
tainty of ±2° in ot through the analysis. 

Second, the accuracy with which wc know the values of the Z 
component ot the dc electric Field is limited bv the performance of 
the PWI Z antenna in these low- to medium-density plasma 
regimes. We have carefully examined the electric field data on a 
case-by-casc basis and selected only those events characterized by 
spin-averaged field profiles which vary smoothly, or with reason- 
ably small amplitude, both immediately prior to and following the 
upwelling event itself, as identified in the ion data. Further, we 
have derived estimates of the magnitude of the uncertainty in the 
values of E z . in the manner described below, and propagated these 
through the analysis. In processing the electric field data, we have 
begun by producing 4-s (spacecraft spin period = 6 s> averages 
/ i J and standard deviations (o 4 s ) about the means of the lb 
-ample s ' data, producing 15 pairs of values per minute. The 
mean values of these two quantities represent l-mtn averages of 
which we have used in the computation of ion bulk 
parameters, and ot a 4 s , which represents a spin period time scale 
error estimate. Further, a 1-min time scale error estimate (cr W) J is 
derived as the standard deviation of the 15 values of E Z>A s about 
their mean. These two error estimates, in addition to a nominal 
uncertainty ot 4 mV m ', have been treated as independent con- 
tributors to the uncertainty in as 

A£ r = \/<i; [ + 5 + (4 ,8) 

4. Results 

4 .1 Results for a Single Pass 

Figure 2 illustrates the nature of the data used in the analysis 
and the results for a single pass on day 60 of 1982. These results 
are not typical, in that this was a particularly intense event, 
observed at unusually low altitude on a day characterized by a sus- 
tained high level of geomagnetic activity. The ion and electric 
Meld data shown in Figure 2 arc plotted against UT. with various 
spacecraft orbital parameters also indicated along the horizontal 
axis. The interval selected shows a DE I pass from the nonhem 
polar cap across the davside auroral zone. 


Figure 2 a shows the RIMS radial head spin time spectrogram 
for O and has been described above in the context of data selec- 
tion. The full time resolution Z component electnc field data and 
the processed Z component data are shown in Figure lb and 
Figure 2c, respectively. The full time resolution electnc field data 
in Figure 2b show not only the steady offset induced by spacecraft 
motion across the magnetic field, but also a distinct modulation 
induced by nonstationary spacecraft sheath structure. Recall that 
since the electric Z antenna is aligned with the spacecraft spin 
axis, spin modulation in this signal is not expected in the presence 
of spacecraft sheath structure unless the structure changes with 
spin phase. Note the large-amplitude, unresolved excursions in 
the electnc field signature near 2327 UT. The long wire electnc 
antenna on DE 1 shows a similar signature, which is quite 
common near the UWI source region. The possible role of these 
large-amplitude and highly structured electric fields in producing 
the upwelling ion flows will be addressed in a future publication. 
The final electric field data product, which is used in computing 
the ion bulk parameters, is shown in Figure 2c. where values of 
A/f. as given in equation (8) are used in determining the size of the 
error bars. In these I -min averages, the contribution induced by 
spacecraft motion has been explicitly removed, and the spin mod- 
ulation removed by the averaging procedure. It is evident from a 
comparison ot Figures 2b and 2c that significant structure, present 
in the unprocessed field data, is lost in the 1-mm data product. 

In Figures 2d through 2/ the bulk parameters (V.,, and n, 
respectively K derived as described above, are presented for each 
°f th e ion species H * , He * , CT . and CT * . Although the field- 
aligned ion velocities and fluxes are generally upward throughout 
the interval, they are most pronounced during the minute extend- 
ing from 2325:30 UT to 2326:30 UT. It is at this time that the 
event is said to maximize in intensity and that the spacecraft is 
thought to be directly downstream of the upwelling ion flow 
source, in the davside polar cleft. 

4.2. Typical Flux . Velocity , and Density 
Levels at Event Maxima 

The ramplike upwelling ton signatures observed as the space- 
craft passes from high latitudes through the davside auroral zone, 
followed by sharp cutoffs at the low-latitude edge, are common 
features ot UWI signatures. These features arc interpreted in 
terms ot a convective effect, whereby the source of the ion up- 
welling exists in a localized zone near the low-latitude edge of the 
event and ions arc convectcd to higher latitudes, dispersing them 
poleward ot the source. From this point of view, the region where 
the upward flux maximizes, near the low-latitude edge of the 
event, may be considered to be most closely associated with the 
source region. The ions observed in this region are characterized 
by the shortest flight paths from source to observer and. as we 
note, arc those with the largest upward field-aliened velocities, 
slower particles having been convcctcd further poleward by the 
time they reach the spacecraft altitude. Wc have therefore taken 
the ion bulk parameters measured at this location to be indicative 
ot the intensity of a given event, with the proviso that there is a 
bias toward larger field-aligned velocities. 

Having selected the analyzed data for 39 events, wc have iso- 
lated the ion bulk parameters at the time in each event when the 
maximum flux was observed. Figure 3 shows a scatter plot of the 
measured maximum flux versus the inferred density for each of 
the tour species: H " . He " . O . and O ^ ~ . Points lying along the 
two lines running through the plot represent upward buik 
velocities ot I km s 1 and 10 km s' 1 . While the measured fluxes 
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! ' o->~ T h ' on data and P WI elect r i c field data are displayed as functions of time for a davside auroral zone pass on dav 60 of 

.h^nj h C °C Sh , 0WS * r ^1° ' spin '" mc counI ra,e spectrogram in standard RIMS format The second and third panels from 
c top show PWI dc electric field data in almost raw and processed formats, respectively (see text) The bottom three panels display 
derived ion species bulk parameters of upward velocity, upward flux, and particle density for several ion species. Note that the O * ion 
densities have been divided by 10 Various orbital parameters, as well as UT. are listed alone the abscissa 


■ind densities vary over several orders of magnitude, the upward 
velocities lend to lie between I and 10 km s 1 for the heavier 
species and near or above 10 km s 1 for the lighter species. Plots 
of the relative occurrence frequency of maximum density, ficld- 
.iliened velocity, and tlux observed in these events for the species 
H .He . O' . and 0" + are presented in Figures 4« through 
4< . respectively In the case of the maximum flux (Figure 4r). the 
flux values plotted along the abscissa have been mapped to a 
common altitude of 1000 km. according to an r ’ scaling law. 
while the densities and velocities are plotted as measured! 
irrespective of spacecraft altitude. This is in contrast to ihe flux 
' .tJucs plotted in Figure 3. which are as measured, irrespective of 


altitude. In Figure 4. the species are differentiated bv the plotting 
symbols used 

Figure 4 a shows that these events are typically oxvgen-rich 
near the source region with typical fractional and absolute O 
ion densities of 0.4 and 1000 tons cm \ respectively. This 
represents somewhat large O' densities for these altitudes 

1.5-2 R f ) as compared with polar cap field lines (see. for 
example, the total electron density results of f\ rsnon ef id, 

I I083|. or the ion species density results of M O. Chandler el al. 
(Observations ol polar ton out t lows, submitted to Journal of 
Ccoffhwmii Rescan h. I9S9). Such large O * densities are indica- 
tive ot increased ionospheric scale heights, probably associated 
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Fig. 3. A scalier ploi ot measured upward ion llux versus measured ion density at the low-laiiiude edee of UWIs is shown for several 

1 . 3and 2 . , The slra.ghf 

show the loci ot points tor ions moving upward at I km s ' (lower line) and 10 km s 1 (upper line) " P 
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t ig. 4. l ine plots showing the frequency of occurrence ot various values 
ot wo ion density, (b) upward bulk velocity, and (c) upward particle tlux. as 
measured at the low -latitude edges tit UWIs. The plotted tlux values are 
normalized to a common altitude of 10* km. 
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unh bulk ionospheric heating ai much lower altitudes. These O' 
density results are consistent with the model polar wind results ot 
Si hunk and Sojka | I9X9|. Our H ' density results are in wide vari- 
ance with those ot Schunk and Sojka. however, with tvpical 
i>hser\ed H densities ot l()0 ions cm at geocentric distances 
between I 5 and 2 R t . as compared with model results m excess ot 
"C\eral thousand ions per cubic centimeter (at 2500 km altitude I. 

The upwelling O' flux (Figure 4< ) is also dominant over that 
associated with the other ions, with maximum sealed flux levels 
which typically approach I0 M ions cm " s ', far in excess ot 
predicted polar wind O r flux levels. On the other hand, the largest 
field-aligned velocities are measured in the H " ions, tor which the 
upward bulk velocity is generally greater than 10 km s' 1 near the 
latitude of the source region. More modest bulk velocities, of the 
order of 3 km s" 1 , are typically observed in the O^ ions. 

We have found that the measured upward field-aligned bulk 
velocity varies systematically with species mass per charge, with 
light ions moving upward with higher velocity than the heavier 
ions. This point is illustrated with Figure 5. where characteristic 
upward bulk velocities arc plotted as a function of species charge 
per unit mass. The distribution of velocities shown in Figure 5 is 
dominated by the FI " ions (at qim - I ) which arc clearly moving 
upward much faster than the other species, although the He" ions 
are aiso significantly faster than the two species (O" and O" 
plotted. The trend is not as strongly supported when comparing 
the O ‘ and O" + ions either in this figure, which is derived from 
the mean values of the bulk velocities observed at maximum flux, 
or on a casc-by-case basis, where the O r is often seen to be 
moving upward with higher velocity than the We have 

shown a power law fit to the data in Figure 5. which yields 

~ j Z ("* km 

v » =l 7 j^( ~ 

where Z is the ion charge state and [mhn p ) is the ion mass, norma- 
lized to the proton mass. Processes yielding energies proportional 
to ion charge but independent of ion mass would give rise to a 
square root dependence of v»| upon < qim). An example of such a 
process is acceleration in a parallel (to B) electric field. However, 
the warm and anisotropic nature of ihesc upwelling plasmas, 
along wuh their high density, suggests that the outflow arises due 
to transverse heating in the topside of an already warm cleft iono- 



l-ig 5. Mean values of ion species upward bulk velocities arc plotted 
versus ihe species charge per unu mass. The velocities plotted arc ihose 
derived irom measurement taken at the lower-latitude edges of LWls. 
w here the upward llux is seen to maximize. The line indicalcs a power law tit 
to die data, with the fit parameters indicated. 


sphere, followed by upward magnetic folding of the ion dis- 
tributions. The analysis of the thermal properties of these events is 
beyond the scope of this report and will be reported in a future 
paper. However, it is to be noted that upwelling ion flows are 
typically warm and anisotropic, with temperatures (transverse to 
B) of several electron volts, as reported by Moore et ai. \ I986/>|. 

Deicourt cr ai. | I989J have modeled the evolution of these ion 
flows using particle rav tracing techniques in realistic gravita- 
tional. electric, and magnetic fields. They have found that the up- 
welling ions generally escape the Earth s immediate vicinity in 
large numbers and travel in an extended plume over the polar cap 
into the geomagnetic tail and plasma sheet, where the heavier 
species are rather dramatically accelerated as they cross the 
neutral sheet. 

4.3. Bulk Parameter Correlations and Dependence 
on Location and Geomagnetic Activity 

The relationship between spatial location and the occurrence of 
ion upwelling has been investigated by Lockwood et ai. |1985|. 
Their study showed that for DE 1 altitudes below 3 R h: . geocen- 
tric. the probability of observing the asymmetric RIMS spin-time 
signature of an upwelling ion event peaks at near 0.7 in the 
prenoon magnetic local time (MLT) sector and at invariant lati- 
tudes (IL) between 75° and 80° (see their Plate 4r). The question 
arises as to whether the intensity of the upward flow shows 
systematic variation with the location of the source region. In 
Figure 6, we present scatter plots of the upwelling flux, meas- 
ured during the minute of maximum intensity within the event and 
normalized to a common altitude of I (XX) km. versus MLT 
(Figure 6</) and IL (Figure 6 b). There is substantial scatter in the 
observations, with several particularly large flux values, and no 
apparent systematic trend. The intensity of the outflow does not 
seem to be correlated with the latitude or MLT at which the obser- 
vations arc made within the source region. 

There is. however, some indication that the ion composition 
within these events varies systematically with magnetic local 
time, with the fractional minor species density increasing toward 
local noon. Figure 7 provides a convenient summary ot the ion 
composition for a number of cases in which all four ion species 
were measured simultaneously. As noted earlier, the bulk of the 
ion plasma consists of O*. comprising about 90# of the ion 
density, while H~ contributes — 10#, and He" and O"" each 
contribute — I #. The upwelling O" flux typically represents a 
smaller percentage (—70#) of the total upward flux, due primari- 
ly to the larger H " bulk velocities. In Figure 7. we show the frac- 
tional ion composition plotted versus the magnitude IMLT - 121 to 
illustrate the existence of a possible systematic relationship 
between proximity to local noon and ion composition within these 
events. Note that values are plotted along the abscissa without 
regard to whether data were collected in the prenoon or postnoon 
sector. The fractional composition of the minor species appears to 
increase with proximity to local magnetic noon. The case for the 
compositional dependence upon MLT is stronger for O " " (linear 
correlation coefficient = -0.7) than for the other species. The 
reason for this dependence is unclear. The O" density (and 
therefore total ion density) appears to decrease with proximity to 
noon, although the data arc scattered (linear correlation coeffi- 
cient |tf|— 0.3), while the absolute minor species densities show 
no significant trend. These apparent trends in the ion composition 
could be anomalous, resulting from correlations between MLT 
and altitude imposed by the DE 1 orbit. Further study of this 
subject will be required. 
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As shown by Lockwood e, al. 1 1985). the invariant latitude c 
the upwellmg ton source region vanes systemat.cally wit 
geomagnetic activity, moving to lower latitudes with incrcasin 
values ot the 3-hour Kp index. The association of this upwellin 
source region with the geomagnetic polar cleft and its svsteman 
latitudinal variation is further illustrated in Figure 8. where th 
invariant latitude at which the maximum upwclling O ~ flux occur 
is plotted versus a measure of the interplanetary magnetic fielt 
i IMF) B. component. The values of B ; plotted were obtained fron 
hourly averages published by the National Space Science Da.i 

Center l NSSDC) World Data Center A for Rockets and Satellites 

hese data are a compilation based on IMP 8 and ISEE 3 observa 
cons. I he plotted values represent weighted land corrected foi 
soiar wind transit timet hourly averages of B. during the houi 


immediately preceding the minute in UT during which the flux 
observation was obtained. For example, in the case of day 71 of 
1982 the maximum O flux was observed at 1934 hours UT. The 
value ot B. associated with that event is taken as 

<»s _ ^l^aivo i-fiXKiil + 26(B ZI mm-moml 

— {Oi 

60 ' ' 

w here S. (jritl(> _ vvn<ll was obtained from the World Data Center data 
base. We have indicated a least squares linear fit to the data, 
which yields a correlation coefficient of 0.72. We have also 
shown (dashed lines) the linear least squares results of Curbary 
and Mem; | I986|. giving the fitted locations of the poleward and 
equatorward cusp boundaries, based on DMSP F-2 and F-4 low- 
energy electron signatures, as functions of B- (measured by the 
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Fig. 8. A scalier plot of the latitude at which the upward O" - flux maxi- 
mizes within UWIs. versus the preceding hourly average of the IMF B. 
component (corrected tor solar wind transit time to the magnetopause). The 
solid line represents a linear least squares fit to the data. The twodashed lines 
show the locations ot the upper- and lower-latitude edges of the clcftlike ( see 
icxt i electron precipitation signatures as functions of B,. These edges arc as 
determined in the study of Carbary and Metu> (1986], 


IMP 8 /GSFC (Goddard Space Flight Center) magnetometer). We 
note here that although Carbary and Meng (1986] use the word 
"cusp," the study they performed applies to a broader region (in 
MLT) which is more properly referred to as the "magnetospheric 
cleft \ Carbary and \feng. 1986, Figure 2 : Newell and Meng, 
1988). The distinction is significant because the cusp, found 
within a region highly localized near noon MLT. as thought to 
connect to the dayside magnetopause, while the cleft, found 
within a broad region which stretches across the dayside, is 
associated with the low-latitude boundary layer. Although, as in 
the Carbary and Meng study, there is significant scatter in the 
data. Figure 8 unambiguously locates the upwelling ton source 
region within the magnetospheric polar cleft, as identified by the 
iow-enorgv electron signature. 

I here does not. however, appear to exist a strong relationship 
between cither the occurrence probability or intensity of up- 
welling ions and the averaged value of the IMF B z during the hour 
immediately preceding the observation of ion upwelline in the 
RIMS data. In Figure 9, we study the occurrence of values of B. as 
observed, both in conjunction with an UW1 observation and 
regardless of UWI observation. All three panels in Figure 9 
feature an abscissa which is composed of B z values binned in I 7 
bins, from - 6 7 to + 67 , with the leftmost bin including all 
observations ot B , less than -67 and the rightmost bin those with 
B z greater than or equal to + 67 . There were a total of 70 
observed UWI events for which the IMF data were available dur- 
ing both the present and previous hours, allowing the weighted 
average described in the previous paragraph to be performed. The 
value plotted along the ordinate ot Figure 9 a represents the 
percentage ot those 70 events which were characterized by a 
weighted hourly average of B z which fell within the given bin. 
The error bar is derived from the square root of the number of 
samples in the bin. For example, 9:r3 of the events were asso- 
ciated w ith a value of B z of less than - 67 . yielding 12.8 ± 4 . 3 % 
ot the 70 samples. In Figure 9 b the occurrence of frequency of B . 
values, lor all available hourly IMF data during the same time 
period as in Figure 9a (days 36-268 of I982L arc plotted for 
comparison. The distribution in Figure 9b is peaked near B , = 0 



IMF B z 

(gamma) 


Fig. 9. The occurrence, during the penod between days 36 and 268 of 
1982. ot various values ot the IMF B. component is studied, boih in relation- 
ship with the occurrence ot UWIs and independent ot such occurrence, uo 
The distribution ot B . values measured during the hour immediately preced- 
ing a RIMS UWI observation. \b) The distribution ol B. for all the available 
data, without regard to RIMS observations, u ) The difference between the 
data in Figures 9 a and 9b These differences should reveal anv significant 
trend in the dependence of the occurrence of dayside ion upwelling upon B. 
All IMF data are corrected for solar wind transit time to the magnetopause 


and falls off quite symmetrically to either side. The large 
occurrences in the extreme bins arc a reflection of the larger bin 
size. Figure 9r displays the differences between the data in the 
above two panels, such that the plotted values are equal to those 
plotted in Figure 9a. less those plotted in Figure 9b. If the occur- 
rence of UWI events were dependent upon B z . we would expect to 
see that reflected most clearly in Figure 9c. If. for example. UWIs 
were more likely to occur for negative values of B z . we would 
expect to see positive differences in Figure 9 c for B z < 0 and neg- 
ative differences for B z > 0 . No such trend is evident, and we 
conclude from these data that UWI occurrence is not dependent 
upon the average value ot B z during the hour immediately preced- 
ing event observation. 

A similar statement can be made regarding the dependence of 
upwelling CT intensity on the averaged B. values, as illustrated 
using the scatter plot ot Figure 10 . Here we have plotted the 
maximum upwelling CT flux, normalized to 1000 km altitude, 
measured during a given event versus the hourly averaged value 
of f° r the h °ur immediately preceding event observation. The 
linear correlation coefficient of -0.14 indicates no significant 
dependence among the parameters plotted. These results are 
somewhat surprising, in view of recent observations of Lockwood 
et al. 1 1988| which suggest that dayside ion upwelling may be a 
direct consequence of magnetic coupling at the subsolar magneto- 
pause. Further study of the relationship of UWI events to 
interplanetary and magnetospheric conditions will be the subject 
of a future publication. 

4 4. Cleft Region Source Strength 

The upward ion flux levels discussed in the previous sections 
are not typical ot the entire dayside auroral zone, but are 
representative ot the most intense portions of the upwelling flows. 
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Fig. 10. Upwelling CT flux is displayed as a scalier plot versus ihe 
immediately preceding hourly averaged value of the IMF B. component, 
corrected for solar wind transit lime to the magnetopause. The flux values 
plotted are the maximum values measured within the UWIs and are normal* 
i/cd to a common altitude of I0 1 km. A linear fit to the data, as well as the 
resulting linear correlation coefficient, is indicated. 


m the immediate vicinity of the source region. We have computed 
bulk parameters from regions downstream (generally poleward) of 
the source region as well. All the measured fluxes, as well as a 
measure of the occurrence frequency of dayside ion upwelling. 
need to be incorporated into any estimate of the total ion outflow 
from this ionospheric source. The expectation value for the up- 


ward particle flux at any location may be expressed as the product 
of the upward flux measured at that location at the time when an 
upwelling ion event is observed and the probabilitv that an up- 
welling ion event will be observed at that location. We have 
binned our ion flux measurements in IL and MIT and combined 
them as a product with similarly binned values of the occurrence 
probabilities reported by Lockwood cr of. | 1985). We have used 5° 
IL bins between 50 : and 9(L and I -hour MLT bins from 0 to 24 
hours. The mean flux in a given bin is taken as 


where Jjj„ represents the mh of a total of A/ flux observations 
(normalized to 1000 km altitude) in the ansular bin associated 
with the indices </./). The expectation value for the flux in the bin 
( /.>> is then given as 

< ^ > »\> “ iJijl X if *ij) ( 11 ) 

where is the occurrence probability reported by Lockwood ct 
al. [1985]. Expectation values computed in this manner arc 
presented in Figure I 1 for the species H T He * . () ~ . and O * * . 
where the value of log M ,t is plotted in polar grey scale 
versus IL and MLT. The expectation values plotted in Fieurc l I 
represent the closest estimate wc can make at this time of the 
mean upwelling source strength tor the various ion species 
shown. The outflowing source region may extend poleward of the 
region shown in Figure I I ; however, such an extension is not 
identifiable within the context of the UW1 signature identification 
on which this work is based. As in the case of the maximum event 
fluxes reported in previous sections. and FI * show the most 
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Fig. 11. The expectation value <•/»> for the upward parallel flux is plotted for the 
grey scale format versus IL and MLT The flux in a given IL-MLT bin is derived as the 
within UWIs and the probability that an UWI will be observed within the bin. 


ion species H ~ . O * . O * . and He ■* . in polar 

product ot the average flux measured in that bin 
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intense outtlows, followed by He* and O * ~ In all eases, the 
source region is rather spatially structureless. 

The fluxes shown in Figure 1 I represent the ion flux expected 
at an altitude of I (XX) km. having been scaled from the observa- 
tion altitudes according to an (nr,,)' dependence, where r and 
are the geocentric radii of the observation point and the l(MH) km 
altitude surface, respectively. The total hemispheric particle out- 
flow rate may be obtained by integrating over the distributions 
shown in Figure II To facilitate the integration, the fluxes arc 
further scaled to the Earth's surface, where 1L and the spherical 
magnetic latitude are equivalent. The total terrestrial outflow- (/) 
from this daysidc source region, for a given ton species, is then 
given as 

f = 2 X<J(/\M j>)> t/a 

ti 

= l d<t> ^ dk <./ l (r u ,A.,4>)> (12) 

where we have approximated the geomagnetic field lines as radial 
and the factor of 2 accounts for source regions in each 
hemisphere. This integration has been numerically carried out 
over the distributions of fluxes shown in Figure 1 1 . The resulting 
outflows are listed in Table I. These numbers may underestimate 
the total outflow, since significant UWI outflows may occur pole- 
ward of regions identified as being within UWls. 


TABLE 1. The Total Rate <J { ,> at Which Ions Flow 
out of the Dayside Cleti Region 


Ion Species 

<i fl >. ions s 1 

FT 

5 x |() :a 

He" 

9 x \{) : - 

O" " 

2 x |() :i 

O" 

2 x I0 :5 


F lows from both hemispheres are included. 


Ihe ion outflow is dominated by O". with 2 x |() : - s ions s 1 
emanating from the dayside cleft region. The H * source strength 
represents —20 7c of the total outflow, while the minor species 
tO" and He") arc at the several percent level or less. These 
results are to be compared with those of Yau et al. (19851 who 
studied major ion outflows at energies such that 10 V ^ Eiq ^ 17 
kV. They reported extensively on outtlows at latitudes including 
the auroral zone and polar cap. distinguishing among flow levels 
on the basis of geomagnetic activity. MLT/IL regime, and solar 
cycle. Many of the results of Yau et a/, j 1985] arc summarized in 
(heir Table I . They found the total net auroral zone ion outflow to 
lie between 2 x I O 25 and 4 x I0 2s ionss 1 forH" and between I 
x 10 ^ and 8 x 10~ s for O depending on geomagnetic activity 
and solar cycle. In our Table 1 . w-c show a total UWI O " outflow 
ot 2 x 1 0- ^ ions s -1 . This is of the same order as that reported by 
Yau et al. during geomagneticallv quiet times. On the other hand, 
our results show a total H " outflow of 5 x 10 :4 ions s" 1 . a factor 
of 5 smaller than that reported by Yau et al. during gcomag- 
neticaily quiet times. We believe these figures apply to ions with 
energies less than or of the order of 10 c\!q. This is based on both 
(he published work of Moore et al, (1986/?) and on other, 
unpublished work, which show typical ion thermal energies 
w ithin these events of 5-10 eV/q. The flux levels in the current 
study are representative of a wide ranee of geomagnetic activity. 
The available data do not permit breaking down the global UWI 


outflow according to geomagnetic activity. It appears, especially 
in the case of the O * . that the UWI source represents a substan- 
tial fraction of the total global ion outflow. 

5. Conclusions 

With the data presented in the sections above we have reported 
the bulk properties of the upwcllmg ion piasma in the davsidc 
cleft region. We have quantified the distribution of ion species 
density, upward bulk velocity, and flux within the cleft source 
region. Additionally, we have combined the bulk flow results 
from a number ot these events with the occurrence probabilities 
obtained by Lockwood et al. { 1985) to give reliable estimates of 
the source strength distribution as a function of MLT and IL. 
Finally, wc have presented a measure of the total cleft region 
source strength in terms of the number of ions per second to flow 
away from the ionosphere. The following specific conclusions 
may be drawn from the work presented above. 

1 . All observed ion species are seen to participate in the ion 
upwclling. Within a given event, the upward flux typically max- 
imizes near the low-latitude edge of the outflow. At that location, 
typical ion species densities, upward field-aligned bulk velocities, 
and upward field-aligned fluxes (mapped to a common altitude of 
1000 km) arc summarized in Table 2. The distributions about 
these typical values have been presented in Figure 4. 


T ABLE 2. Typical Ion Bulk Parameters as Measured Within, but Near the 
Lower-Latitude Edge of Upwclling Ion Events 


Species Bulk Parameter 

O* 

H" 

He" 

()" " 

n. cm ' 

600 

50 

15 

7 

i’ll, km s 1 

3 

16 

7 

4 

7 f |,* particles cm ' s 1 

6 x 10 * 

2 x 10* 

I0 7 

6 x I0 (> 


^Scaled to a common altitude of 1000 km. 


2. There is a systematic dependence of the upward field- 
aligned ion bulk velocity upon the species charge-to-mass ratio, 
with the velocity increasing approximately with the square root of 
<■( w- This dependence is illustrated in Figure 5. 

V The intensity ot UWIs. as measured by the maximum up- 
ward O" flux observed within an event, docs not show significant 
systematic variation as a function of MLT or IL. This is indicative 
of a rather uniform outflow from the dayside cleft reeion. This 
point is illustrated in Figure 6, where scatter plots of the max- 
imum normalized O" flux versus MLT and IL are presented. The 
point is further illustrated, with respect to all measured flux values 
for the species H * , He". O". and O"*. in Figure II. where 
the expectation value of the normalized upward flux is shown 
plotted in polar form. 

4. The ion density composition of the outflow shows a weak 
dependence upon MLT. however, with minor species, particular- 
ly 0~ representing a larger fractional composition near magnetic 
local noon. This point is illustrated in Figure 7. 

5. The statistical association of UWIs with the dayside mag- 
nctosphenc cleft, as identified by low-energy electron signatures 
| Carbary and Mem* , 1986). is illustrated in -Figure 8. This associ- 
ation points up the likelihood of a relationship between the 
phenomenon of dayside ion upwelling and the solar wind- 
magnetosphere interaction. 

6 There is. however, no clear-cut relationship between cither 
the occurrence or intensity of UWIs and the immediately preced- 
ing hourly averaged value of the IMF B z component. These facts 
are illustrated in Figures 9 and 10. respectively. As the simplest 
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measure ot the degree of electromagnetic coupline between the 
magnetosphere and solar wind, one might expect to'find increased 
occurrence or intensity ot UWIs with increasingly negative 8. ,1'. 
lor example. UWIs were directly associated with variable aspects 
Ot t c large-scale auroral convection morpholouv or with the 
occurrence ot magnetopause flux transfer events.' both ot which 

depend upon «. The results presented here indicate that neither ,s 
the case. 

Our estimate of the ion outflow rates (including both 
hemispheres: trom this davside cleft source have been presented 
m Fable 1 These estimates arc given based on data collected 
between days 36 and 137 of 1982. and almost all were from 
northern hemisphere passes. Thus, they are generally representa- 
tive ot spring conditions during an active portion of the solar 
cycle. Virtually all levels of geomagnetic activity are represented, 
n the case of C> + and H". these numbers validate the cleft ion 
fountain contributions used by Chappell et al. [1987] and 
Delcoun et al. [ 1989) in demonstrating the significance of the 
ionosphere as a magnetospheric plasma source. 
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Statistical survey of pitch angle distributions in core (0 
50 eV) ions from Dynamics Explorer 1: Outflow in the 
auroral zone, polar cap, and cusp 


B. L. Giles , 1 C. R. Chappell , 2 T. E. Moore , 1 R. H. Comfort , 3 and J. H. Waite Jr . 4 

Abstract. Core (0-50 eV) ion pitch angle measurements from the retarding ion mass spec- 
rorncter on Dynam.es Explorer I are examined with respect to magnetic disturbance, invariant 
luti ude magnetic local time, and altitude for ions H*, He*, O*, M/Z=2 (D* or He**) and O** 
Included are outflow events in the auroral zone, polar cap, and cusp, separated into altitude re- 
c ions below and above 3 R E . In addition to the customary division into beam, conic and up- 
inn the h.gh-latitude observations fall into three categories corresponding to 

This b s!nW dS h | 3t afe u f S ? an ’ (2) com P arable t0 - or (3) faster than that of the spacecraft 
Jnhpri P a ° ng WUh the a t,tude part,non ’ serves to identify conditions under which iono- 

spheric source ions are gravitationally bound and when they are more energetic and able to escape 
<o the ou ter magnetosphere. Features of the clef, ion fountain inferred from single event stuST 
are clearly identifiable in the statistical results. In addition, it is found that the dayside pre-noon 

level aniUhP T h 50 f fCe ° C ,T Vel0C,ty low - ener gy '°ns regardless of species or activity 

to vhv d Th d yS ' de . ?° 0n ^ ° r aUr ° rai ZOne ‘ h^ 0 " 165 an additional source for increased 
activity. The auroral oval as a whole appears to be a steady source of escape velocity H* a 

steady source of escape veloc.ty He* ions for the dusk sector, and a source of escape velocity 
heavy ions for dusk local times primarily during increased activity. The polar cap above the au- 

neart n°h e “ 3 f° nS,Stent . source of lo w-energy ions, although only the lighter mass particles ap- 
pear to have sufficient velocity, on average, to escape to higher altitudes. The observations sup- 
port two concepts for outflow: ( I ) The cleft ion fountain consists of ionospheric plasma of I -20 

Lr rgy Stre r ng T ard ,nt0 the magnetos P her e where high-latitude convection electric 
fields cause poleward dispersion. (2) The auroral ion fountain involves field-aligned beams 

which flow out along auroral latitude field lines; and, in addition, for late afternoon local times 

£ ! P CC add ‘ t,onal acceleration such that the ion energy distribution tends to exceed the 
detection range of the instrument (>50-60 eV). cauccu me 


1. Introduction 

At polar latitudes, because the Earth's dipole field is nearly 
vertical, charged particles in the ionosphere created by pho- 
loionization or particle impact ionization are relatively free to 
move upward along the magnetic field lines and populate 
higher-altitude regions of the magnetosphere. Since the dis- 
covery of energetic O* ions by Shelley et at. [1972], ion com- 
position measurements have revealed that a significant 
amount of plasma in the Earth's magnetosphere is of iono- 
spheric origin [e.g., Chappell et at., 1987; Moore et at., 
1989; Moore , 1991], The importance of this terrestrial source 
is confirmed by the presence of energized ionospheric ions 
throughout the magnetosphere, such as the ring current and 
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piasma sheet [Johnson et at. , 1977; Geiss et at., 1978; 
Lennartsson et at. , 1979; Batsiger et at., 1980; Lundin et at.; 
1980, Peterson et at., 1981; Sharp et at.. 1982; Eastman et 
ai, 1984; Lennartsson and Shetley , 1986], the magnetopause 
and magnetosheath region [Peterson et at. , 1982; Hultqvist , 
1983], the magnetotail lobes [Sharp et at., 1981], and the 
plasma mantle and magnetotail boundary layer [Frank et at., 
1977; Hardy et at., 1977; Lundin et ai, 1982a. b; Candidi et 
at., 1982]. In addition, it has been shown that a significant 
portion of these ionospheric-source ions have energies below 
1 keV [Yau et at. , 1934; Giles et at. , 1988; Delcourt et at., 
1988]. These lower-energy particles are very sensitive to 
small fluctuations in local magnetic and electric fields. They 
interact with charge carriers for magnetospheric and iono- 
spheric currents and operate as a permeating medium in which 
plasma waves are generated and propagated to other parts of 
the magnetosphere thereby influencing larger-scale interac- 
tions. Thus information about the dynamics of this plasma 
component is an effective tracer for broader scale magneto- 
spheric processes. 

In understanding the transport and dynamics of charged par- 
ticles, information about the particle pitch angle distribution 
(PAD), the arrangement of directional flux with respect the lo- 
cal magnetic field direction, is an invaluable aid. For exam- 
ple, case studies by Sharp et ai [1977, 1979] exploited Held- 
aligned pitch angle distributions to infer the presence, and 
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vertical extent of, parallel potential drops within an auroral 
acceleration region. Also, occurrence statistics of distinct 
PAD classifications have been used to examine possible 
plasma source locations, acceleration regions, and the subse- 
quent evolution of plasma distributions during transport [e.g., 
Nagai et al., 1983; Yau et al ., 1984], The lowest-energy ions, 
in particular, have been shown to exhibit complex pitch angle 
characteristics such as lield-aligned beams, conics, and 
trapped (pancake) distributions [ Chappell , 1982; Chappell et 
al. , 1982a]. 

The retarding ion mass spectrometer (RIMS) experiment on 
Dynamics Explorer 1 (DE 1) represents a major effort to thor- 
oughly characterize the lowest-energy plasma component by 
measuring ions ranging from 0 eV (or spacecraft potential) up 
to 50 eV throughout almost the entire volume of the inner 
magnetosphere, from the topside ionosphere out to 4.67 R E . 
The RIMS effort has, in its lifetime, revealed previously unob- 
served ionospheric plasma sources to the magnetosphere and 
called into question accepted theories concerning the relative 
contributions of the internal and external sources to magneto- 
spheric plasma content. Specific to the RIMS instrument are 
observations of the supersonic polar wind [Nagai et al., 
1984], a source of low-energy heavy ions from the polar cap 
and the cleft region [Waite et at'\ 1985; Lockwood et al. , 
1985a, b], upward fluxes of N*, N** [Chappell et al., 1982b] 
and of molecular ions N 2 *, NO*. and 0 2 * [Craven et al., 1985]. 
It is this data set that will be used to obtain pitch angle distri- 
butions of the individual charged particle species. 

One of the most intriguing topics in magnetospheric 
physics research that will be addressed is the question 
regarding the source of particles populating the inner magne- 
tosphere. Concentrations and relative dynamic behaviors of 
H*. He**, and O* are commonly used as indicators of a solar 
wind or ionospheric source. The RIMS mass spectrometer is 
capable of resolving several mass peaks, including H\ He*, 
OL O**, and one at a charge-to-mass ratio of 2. Because RIMS 
provides no direct information on charge state, these latter 
measurements could conceivably be H 2 + , He**, or D*. Breig et 
al. [1987] and Breig and Hanson [1991] have reported on deu- 
terium ion Hows in the thermosphere. In this case the instru- 
ment is a predecessor of RIMS which again, yields the M/Z-2 
mass peak with no charge information, but additional correla- 
tive evidence with H* determines these samples to be D*. 
Young et al. [1977, 1982] and Comfort et al. [1988] report 
observations in the relatively cold outer plasmasphere that are 
most likely terrestrial source He**. In the present work, be- 
cause the DE 1 orbit samples low altitudes where D* measure- 
ments are likely and samples higher altitudes where low- 
energy He** measurements are possible, the RIMS 
measurements will be referred to as being of the M/Z-2 ion. It 
is thought, because of the RIMS energy range, that any He** 
measurements would represent doubly-ionized helium produced 
in the ionosphere or plasmasphere, not cooled solar wind 
He**. Regardless though, of whether the M/Z-2 measurements 
are representative of He** or D*. the key issue is whether 
evidence exists for these ions, and of H*, He*, O*. and O**, to 
have an ionospheric source and then flow outward, such that 
they populate both inner and outer portions of the 
magnetosphere. 

In section 2 we describe the RIMS instrument and the DE 
orbit, and outline the pattern recognition technique used. 
Section 3 describes the classification of the data in terms of 
pitch angle distributions and details the consequences of ob- 
serving low-energy particles with a moving spacecraft. 


Section 4 contains occurrence probability data with respect to 
magnetic disturbance, invariant latitude, magnetic local time 
(MLT), and geocentric altitude for invariant latitudes above 
60°. Section 5 discusses these data in relation to the transport 
and dynamics of plasma in the Earth s inner magnetosphere. 
Section 6 is a summary. A paper (Giles et al.. Statistical sur- 
vey of pitch angle distributions in core (0-50 eV) ions from 
Dynamics Explorer 1: Map ot the magnetosphere’s low- 
energy ion component, unpublished manuscript, 1994) 
presents statistics for invariant latitudes below 60°. 
examining trends with disturbance, magnetic latitude. MLT 
and L shell. 

2. Instrumentation and Data Analysis 
Considerations 

The survey covers the ions routinely measured by the DE 
1/RIMS instrument (H*, He*. O*, M/ Z=2, and O**) in the energy 
range 0-50 eV. The observed count rates for each ion were 
classified according to their pitch angle distribution as deter- 
mined from the radial sensor head of the instrument as it ro- 
tated through a full range ot pitch angles in each 6-s satellite 
spin period. As the spacecraft spins, the radial head (mounted 
perpendicular to the spin axis) responds to ion (luxes for 
nearly the full range of angles (-180° to +180°) with respect to 
the spacecraft ram direction. With knowledge of the Earth s 
magnetic field in relation to spacecraft ram. the RIMS spin 
angle* response corresponds closely to a pitch angle distribu- 
tion (detailed in section 3). The instrument and its operational 
characteristics are described by Chappell et al. [1981], Fields 
et al. [1982], and Olsen et al. [1985]. 

The elliptical polar orbit has perigee of 675-km altitude and 
apogee of 24,875-km altitude. The motion of the Earth around 
the Sun causes the orbit plane to drift westward at a rate of I 
hour MLT every 15.4 days (covering 24 hours in 12 months); 
the orbit line of apsides drifts 0.328° geographic latitude each 
day (moving from one pole to the other in 18 months). All 
available data, from October 19. 1981, to December 31, 1984 
(a full precession/rotation cycle), are used to provide nearly 
complete local time and latitudinal coverage out to 4.67 R E . 
However, the 3 to I ratio between the drift of the line of ap- 
sides and the westward drift of the orbit plane does result in 
specific areas of poor sampling in local time for a given alti- 
tude. For example, in the 0 to 6 MLT sector, midaltitudes in 
the auroral latitude regions are not well sampled. In addition 
to regions not reached by the orbit, the instrument is not al- 
ways operated down to the DE 1 perigee of 675 km. At low al- 
titudes the plasma density can be quite high and the instrument 
often enters a shutoff mode to protect the particle-detecting 
channeltrons from excessive count rates. 

The instrument operational modes allowed programmed 
sampling of the mass range such that the major masses were 
sampled more frequently and constitute a larger portion of the 
database. Table 1 lists the number of 1-min spin angle distri- 
bution samples for each species above 60° invariant latitude 
(ILAT) in four MLT quadrants and two geocentric distance lev- 
els. Over 500.000 samples were obtained for the H\ He*, and 
O* ions and over 400,000 samples for the M/Z= 2 and O** ions. 
Because ot the varying sample rate, the statistical calculations 
are carried out for each ion species separately; the probabili- 
ties presented are characteristic of the population of individual 
ion species, not ot the ion population as a whole. Fifty-eight 
percent ot the data were obtained during periods of relatively 
low activity (Kp < 3) and 42% during high activity (Kp > 3). 
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It should be noted that, for the 0 to 50 eV energy range, 
floating spacecraft potentials can prevent the detection of the 
I o west-energy ions. The RIMS team dealt with this spacecraft 
charging by the addition of a biasable voltage plate at the en- 
trance aperture programmable to 0, -2 -4, or -8 V relative to 
the spacecraft potential. The present data set includes all 
available data, regardless of aperture bias setting. It is possi- 
ble that the data set does not include particles with ram energy- 
per-charge less than the potential of the spacecraft. 

Por each ion the raw RIMS data are read and the counts-per- 
accumulation period (12 ms) are averaged over the time resolu- 
tion period for each spin angle bin. For this study, I -min time 
averages with 5° spin angle bins were chosen. This 5° spin 
angle bin refers to the center of the view cone and does not 
imply 5° resolution in the measurements. A subroutine scans 
the counts-versus-spin-phase-angle array to identify and char- 
acterize. in terms of centroids, widths, and skews, local max- 
ima in the distribution with respect to the magnetic field direc- 
tion. Next, a filtering pattern recognition routine produces 
the spin angle distribution description for each minute of 
RIMS data. When combined with orbit parameters and mag- 
netic activity indices, this is the source for statistical analy- 
sis. Further description of the data analysis may be found in 
the work by Giles [1993], 

The relations used to calculate the centroid and coefficient 
ot skewness, taken from Sever (1984], apply to discrete data 
grouped into a distribution having spin angle bin classifica- 
tions 0,, I = 1, 2, .... k. and corresponding measurement val- 
ues (in this case, instrument count rate per spin angle bin), 

i=l,2 k. The centroid is an arithmetic mean given 
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where W 2 and M } are the second and third moments about the 

mean of the sample. The rth moment about the mean M is 
given by r ’ 

<«> 

p.T l “ sir * ca ‘ i " n °f Spin Angle Data in Terms of 
Pitch Angle Distributions 

3.1. Relationship to Pitch Angie Distributions 

Strictly speaking, spacecraft spin angle distributions ob- 
served in situ are not true pitch angle distributions as custom- 
arily defined. The measured spin angle is defined in the refer- 
ence frame of the spacecraft, which is different from the 
plasma rest frame due to spacecraft motion (2 to 8 km s 1 de- 
pending on spacecraft altitude). This difference in reference 
lmle consec l uence provided the thermal speed. 

( kl/M) , of the plasma is large compared to the spacecraft 
ram speed. For an average spacecraft speed of 5 km s 1 a fac- 
tor-of-5 ratio of thermal speed to spacecraft speed requires 3 5 
eV for H*. 1 3 eV for He*, and 52 eV for O*. 

Because RIMS responds to ions with energies considerably 
below this threshold, the RIMS spin angle data contain infor- 
mation about both the angular distribution of the plasma in its 
own frame of reference (the conventional pitch angle distribu- 
tion) and the drift motion of the plasma relative to the space- 
craft. These two aspects of instrument response tend to domi- 
nate independently - one at times obscuring the other - de- 
pending on the plasma's thermal speed and its relation to the 
spacecraft and plasma drift speed. Figure I illustrates the ex- 
tremes of this motional influence. 

Figure I shows simulated RIMS spin angle responses to a 
bi-Maxwell.an plasma distribution ( TJT „ = 5) for four combi- 
nations ol plasma thermal and drift speeds. The simulation 
adopts a numerical quadrature approach to the solution of the 
response integral, using a sum over aperture area, solid angle, 
and energy and weighted with the appropriate instrument re- 
sponse factor to obtain the correct energy-angle dependent re- 
lationship between detector count rate and directional flux. 
The parameters illustrating the limits of response were chosen 
for ease of understanding but without strict adherence to phys- 
ical likehood in any given region of the polar ionosphere. In 
each case the spacecraft speed is set at 5 km s' 1 perpendicular 
to B. Shown beside each spin angle plot are two phase-space 
schematic diagrams to clarify the orientation of the various 
spacecraft and plasma components. Both display velocity 
space, one in the spacecraft reference frame, the other in the 
frame of reference of the plasma. 

In Figure la the H* ions have the perpendicular thermal 
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3 RIMS Modelled Response: V th > V sc 



R * Response direction 

RIMS Modelled Response: V th < 



R * Response direction 


Figure 1. Simulated retarding ion mass spectrometer (RIMS) spin angle responses to bi-Maxwellian 
plasma distributions (T 5). Also shown are phase space diagrams, drawn schematically, locating the 
various spacecraft and plasma components. Note that instrument response is to ions with velocities anti- 
parallel to the detector look direction, (a) H + ions with perpendicular thermal speed component greater than 
spacecratt speed by a factor of 5. The top panel assumes a stationary plasma and the bottom panel assumes 
Vq = 25 km s l , parallel to B. (b) 0 + ions with perpendicular thermal speed a factor of 5 below the spacecraft 
speed and with same drift speed assumptions as Figure la. For plasmas with thermal speeds smaller than the 
spacecraft and drift speeds, relative drift motion dominates the response. For plasmas with thermal speeds 
larger than the spacecraft and drift speeds, the true plasma frame distribution function emerges. 


speed component greater than the spacecraft speed by a factor 
of 5 (Vth - 25 km s' 1 and V sc = 5 km s 1 ). In the top panel the 
RIMS response is to a stationary plasma (drift speed V D = 0); 
the angular characteristics of the distribution are clearly visi- 
ble as peaks in the spin angle directions perpendicular to the 
magnetic field (0° and 180°). The relative motion of the 
spacecratt (5 km s 1 with respect to the stationary plasma) ap- 
pears as an elevation in count rate in the direction of space- 
craft motion (0°). The bottom panel assumes V c = 25 km s* 1 , 
parallel to B. In this case, bulk drift speed comparable to 
thermal speed, the angular characteristics of the distribution 
have become nearly unresolvable although the additional drift 
speed of the distribution may be identified by the shift from 


ram in the direction corresponding to an addition of the space- 
craft and drift velocity vectors. (Without regard to the rela- 
tionship of V ( h to V sc , simulation of the instrument response 
shows that angular features are resolvable for TJT\\ ratios 
greater than -1.5). 

In Figure lb the response is for a bi-Maxwellian distribu- 
tion of O* ions with perpendicular thermal speed a factor of 5 
below the spacecraft speed ( V, h ~1 km s' 1 and V sc = 5 km s* 1 )- 
In the upper panel, even though V D = 0, the angular character- 
istics of the bi-Maxwellian are not resolvable within the de- 
tection thresholds of the instrument and the distribution ap- 
pears as a single peak at the spacecraft ram direction. The bot- 
tom panel assumes V D = 5 km s 1 parallel to B, the additional 
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drift speed of the distribution again apparent as a shift fro 
ram. 

In summary then, for plasma with thermal speeds mui 
smaller than those of the spacecraft and drift speeds, the rel 
live drift motion dominates the response, leading to 
rammed flux peak in the direction of relative wind as seen 
the spacecraft frame. In this case the angular information 
given by subtle features of the detailed shape of the ramme 
peak but is often indeterminate. For plasma with thermi 
speeds much larger than those of the spacecraft and dri 
speeds, the true plasma frame distribution function emerge 
clearly, although the relative drift motion, appearing as sma 
asymmetries in the angular information, may be difficult to re 
solve. The RIMS data included in this survey span this ent.r 
range ol extremes, although the specific set of plasma parame 
ters contributing is mass dependent. 

Figure 2 illustrates this for H* and O* by plotting plasma pa 
rameters covered by the data set as a function of ion Macl 
number. Here ion Mach number is defined as the ratio of th. 
distribution drift velocity (in the Earth frame) to the ion ther 
mal speed (2 kT/M)' n . Roughly speaking then, it is for sub- 
sonic (thermal speed greater than drift speed) ion Mach num- 
bers that the angular characteristics of the distribution will be 
resolvable. Figure 2 shows that for drift velocities less than 
5 km s ' H *' wn™* are subsonic for thermal components as 
low as 1.0 eV. Hence the angular characteristics of the H* 
samples in the data set should be easily resolvable except for 
the coldest plasma with strong (V D > 15 km s' 1 ) drift veloci- 
ty. O- outflows are subsonic at 15 km s ' for thermal compo- 
nents above 10 eV and at 8 km s ' for thermal components 
above .1 eV. Therefore, for heavy ion outflows with tempera- 
tures below 5-10 eV. the angular information will frequently 
ne indeterminate. 

Similarly, the drift motion effect is resolvable in the super- 
some Mach number range. Because for small thermal compo- 


nents it is the relative difference between drift and spacecraft 
velocity that is important. Figure 2 can be reinterpreted in 
terms ol a zero ion drift velocity. The shaded region then 
serves to delineate the range of spacecraft velocities with alti- 
tude (-7.6 km s ' at 1.5 R E geocentric, -4.2 km s' ' at 3 % and 

~. " km s at 4 5 r e> an d therefore the ion Mach range for a 
discernible relative drift motion effect to the spin angle mea- 
suremems (i.e.. not dominated by the angular response). In 
the shaded region the relative difference between drift and 
spacecraft velocity is small, so that light ion (H" ) outflow 
measurements with thermal energies below -0.4 eV will be 
dominated by the relative drift motion - a range for which de- 
tection is difficult due to positive spacecraft potentials For 
heavy ions (O') the drift motion effect is significant for ther- 
mal energies up to 5 eV. and so will be a frequent spin angle 
feature in the heavy ion outflow samples 


3.2. The Classification Scheme 

Our approach is to separate the spin angle data into cate- 
gories representing distinct physical behaviors, both for an- 
gular characteristics (the conventional pitch angle distribu- 
tion) and on the basis of apparent drift speed relative to the 
spacecraft. As a consequence to this, since the local ion es- 
cape speed is near to that of the spacecraft, the apparent drift 
speed allows us to categorize plasma behaviors loosely in re- 
lafion to the local escape speed. 

Figure 3 shows data measured by RIMS, and they are given 
as examples of the distribution types featured in this paper. 

I he data are plotted as accumulated counts per sample versus 
spin phase angle. In addition, the detector pitch angle is plot- 
ted as a (unction of spin angle. The counts per sample is pro- 
portional to the integral ion flux for the energy range 0 to 50 
e with approximate conversion ratios of counts to integral 
flux of 3.1xlO J counts cm- s 1 ster 1 for H*. I.9xl0 4 counts 
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DE*1 RIMS 81/287 20:37 - 20:38 DE-1 RIMS 82/071 19:34 ■ 19:35 



Spin Angie (Deg) 


Spin Angle (Deg) 


DE-1 RIMS 81/295 13:38 « 13:39 DE-1 RIMS 81/285 22:50 • 22:52 



Figure 3. Examples of DE 1 RIMS spin angle distributions typical of high-latitude outflow. Approximate 
conversion factors from count rate to integral ion flux are given in the text. Illustrated are: (a) unidirectional 
field-aligned beams, (b) uni-directional field-aligned conics, (c) asymmetric distributions characteristic of 
the upwelling ion signature, and (d) rammed distributions, which illustrate the effect of comparable space- 
craft and plasma outflow velocities as discussed in section 3. 


cm' 2 s' 1 ster' 1 for He\ and 4.4xl0 4 counts cm 2 s l ster 1 for 0 + . 
These factors, which reflect the mass-dependent angular re- 
sponse, are based on preflight and in-flight calibration proce- 
dures; the ratios are time-dependent and are offered here only as 
a general approximation. 

The unidirectional field-aligned beam distributions (Figure 
3a) are characterized by a single flux peak certered at the field 
direction (here near 180° pitch angle). These distributions 
may be shifted from the field direction due to a combination of 
convection, drift, and spacecraft velocity effects. The impor- 
tant distinction is that these distributions are ordered by pitch 
angle (magnetic field direction) instead of by spin angle (the 
velocity vector of the spacecraft) and therefore correspond to 
field-aligned flows with velocity greater than that of the 
spacecraft. The rammed distribution (Figure 3d) is most com- 
monly identified in the data by a dome-shaped symmetric peak 
centered at or near the satellite ram direction. From the argu- 
ments above, the rammed distributions correspond to plasma 
flows that are slower than the spacecraft. The unidirectional 
field-aligned conic (Figure 3b) are those having flux peaks be- 
tween the directions perpendicular and parallel to the magnetic 
field direction. In this example, the spin angle distribution 
appears as a large single peak but with a distinct drop in count 
rate at the field direction (180° pitch angle). In other exam- 
ples, the conic can appear as two separate peaks. The asym- 
metric distributions (Figure 3c) appear as a single peak in 
count rate near the spacecraft ram direction but with a notable 
asymmetry toward the upward magnetic field direction. RIMS 
samples that have a skew coefficient (equation (5)) exceeding 
0.75 toward the upward magnetic field direction will be catego- 


rized asymmetric. In the example shown the coefficient of 
skew is 0.798 defined about the centroid position (14.6° in 
this example) for the portion of the distribution (61.7° in spin 
angle width) above the count level equal to l/<? 2 times the max- 
imum. It is this data signature that is associated with the up- 
welling ion events first identified by Lockwood et ai [1985a). 
Distributions fall into the intermediate category when a single 
peak in the count rate falls between the field and ram direc- 
tions. For this reason, they cannot be readily assigned to ei- 
ther category and tend to be the samples in which the plasma 
and spacecraft velocities are more nearly equal. Within the 
RIMS data set, downflowing field-aligned beam distributions 
primarily consist of ions falling under the influence of grav- 
ity, the velocity of which is expected to be only a few kilome- 
ters per second. The unidirectional field-aligned downflows 
category is a subset of the rammed and intermediate categories 
and contains the ion downflows for which a clear downward 
drift velocity is resolvable. 

The remaining distribution classifications covered in this 
paper, nothing observed and low counts, are to classify data 
signatures that do not fall into one of the categories defined 
above. The data intervals that have counts below the 2 counts- 
per-sample noise level are assigned to the ’nothing observed" 
category. The data intervals that have counts above the noise 
level, but without identifiable maxima are assigned to the 
"low-count" category. A later paper will examine distribution 
types primarily identified with midlatitudes (i.e., counter- 
streaming and trapped); all distribution types are included in 
the counting statistics so that the porportions reported are 
representative of the full population. 
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Figure 4. Est. mated bulk drift velocities for samples in the rammed, intermediate, and beam PAD cate 
ton'nf fh 3 Un ?°T h ° f , SpaCeCraU Ve,0city (and consequently, as a function of altitude which is noted at the 
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de f v wh ch PAD "* ™» "I*— by observed spin angle distribution serves to 

“ "f A ° cate 8° nes conta,n gravitationally bound ions and which contain ions that are more en- 
ergetic and thus able to escape to outer magnetospheric regions. 


Figure 4 quantifies the inferred outflow drift velocities asso- 
ciated with three of the above categories by providing a rough 
estimate of the bulk drift velocity for each sample as a func- 
tion of spacecraft velocity (and consequently, as a function of 
altitude). The heavy black line marks the ion drift velocity 
necessary to escape Earth's gravitational potential. The inten- 
tion is to bracket the range of plasma drift speed that result in 
a given sample being assigned to one category or another. 
Other studies by DE investigators [Chandler et ai , 1991; 
Pollock et al. % 1990] report detailed velocity statistics for po- 
lar wind and cleft ion fountain outflows. 

The calculation assumes the plasma drift velocity to be 
aligned with the magnetic field (i.e., no convection) and that 
only the spacecraft velocity component transverse to the field 
may bias the observed distribution toward the direction of 
spacecraft motion. We use the approximation that the mag- 
netic field lies in the spacecraft spin plane . This is generally 
valid within 12° for the present data set (and within 5° for 91% 
ot the samples). Provided the spin angle corresponding to the 
centroid of the observed plasma distribution (a) and the spin 
angle corresponding to the magnetic field direction (p) do not 
coincide, the plasma drift velocity may be calculated from 

x/ _ ^spacecraft s ' n (°0 

piasma- sin (/3 - or) (7) 

Figure 5 is a sketch of the vector measurements and angles 
relating to this calculation. For the rammed category all out- 
flow velocities, with few exceptions (all H + ), fall below the re- 
quired plasma escape velocity across the entire altitude range; 
consequently, the samples represented by this category, with- 
out additional energization, must either reenter the ionosphere 
or he trapped by magnetic mirror effects. The intermediate 
category is a mixture of escaping and nonescaping H + sam- 


ples, and contains mainly nonescaping OT The beam cate- 
gory contains the majority of escape velocity samples for 
both ion species, although there is a significant portion of 
relatively fast, 1 but nevertheless nonescaping, 0 + at the 
highest altitude (lowest spacecraft velocity) bins. 

At high latitudes, then, the field-aligned outflows fall into 
three categories (rammed, intermediate, and field-aligned 
beams), which, for the energy range detectable by RIMS, cor- 
respond to ion flow speeds that are (1) less than/(2) compara- 
ble to, or (3) much faster than that of the spacecraft. This sep- 
aration by observed spin angle distribution serves to identify 



Figure 5. Sketch of the vector measurements and angles for 
the bulk drift velocity calculation shown in Figure 4. The 
value a is the centroid of the observed plasma distribution, 
and P is the spin angle corresponding to the magnetic field di- 
rection. The calculation assumes the plasma drift velocity to 
be aligned with the magnetic field (i.e., no convection) and 
that only the spacecraft velocity component transverse to the 
field may bias the observed distribution toward the direction of 
spacecraft motion. 
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Figure 6. (a) The percentage ot samples, during which RIMS was programmed to accept a given ion, the 
instrument was able to detect that ion. Shown are seven Kp bins for H*, He*, 0 + , M/ Z=2, and 0**\ The occur- 
rences are calculated for invariant latitudes above 60° without regard to PAD category, altitude, or local time, 
(b) In six panels, occurrence with altitude for asymmetric, unidirectional field-aligned conics, unidirectional 
field-aligned beams, rammed, intermediate, and nothing observed; the calculations include only data above 

60 ILAT. The prominent feature to note is the separation of the various distribution types into distinct alti- 
tude ranges. 


which ions will be gravitationally bound and which are more 
energetic and thus able to escape to outer magnetospheric re- 
gions. 

4. Statistical Results 

4.1. Occurrences With Kp and Altitude 

The RIMS instrument, when programmed to detect a selected 
mass species, may or may not detect appreciable ion counts 
above the noise threshold within the 0-50 eV energy range. 
Figure 6a summarizes, as a percentage of the total number of 
samples, when RIMS measured ion counts above a 2 counts- 
per-sampie threshold. Data are shown lor seven magnetic ac- 
tivity level (Kp) bins for the ions, H*. He*, O*, M/Z= 2, and 
0 ++ . The occurrences are calculated for invariant latitudes 
above 60°, without regard to altitude, local time, or the cate- 
gories of pitch angle distributions defined above. The figure 
simply indicates in what percentage of samples, during which 


the instrument was programmed to accept a given ion, the in- 
strument was able to detect that ion. 

Counts of 0-50 eV H + are measured above the detection 
threshold in 75-78% of the samples, with no significant varia- 
tion with magnetic activity. In contrast, occurrences of the O* 
ions are strongly dependent on the activity level. At the low- 
est Kp bin, O* ions are measured in 33% of the samples; this 
percentage rises steadily with activity, exceeding the occur- 
rences of H + for the highest Kp bin. The He* ion occurrence 
probabilities lie between those of H* and O* and point to a 
weak dependence on activity. Displaying a pattern similar to 
that of O*, O** appears in about 19% of the samples for low Kp 
levels and increases to 41% during times of high Kp. M/Z= 2 
ions are seen in 20-29% of the samples, with no clear trend 
with activity. 

Figure 6b has six panels showing the percentage of occur- 
rences as a function of altitude for: (1) asymmetric distribu- 
tions. (2) unidirectional field-aligned conics, (3) unidirec- 
tional field-aligned beams, (4) rammed, (5) intermediate, and 
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ItT? ° b , SerVed; ‘ he calculations delude only data above 
ou ILAT and are sorted without regard to local time. The 
probability calculations are carried out for each ion species 
separately. The prominent feature to note in the altitude pro- 
files is the separation of distribution types into distinct alti- 
tude ranges. The asymmetric, intermediate, and rammed distri- 
butions appear with greatest frequency at altitudes helow 2.25 
^ £ °r a |' ,ons - Unidirectional conics, most clearly identified 
^ f^ ples ' peak in Probability at midaltitudes between 
1/5 and 3 R E . Unidirectional beams tend to populate the 
higher end of the altitude range; they reach maximum fre- 
quency above 2.5 R E for H- and He- and 3.5 R E for O- The 
nothing observed category, samples with counts below the 
no.se level, also peaks in probability at higher altitudes for 
He , O . and 0~. Statistics for the M/Z=2 PAD types, not eas- 
i y viewed in this format, display relatively higher frequencies 
at the lowest altitudes with the exception of intermediate 
distributions which occur with uniform frequency across the 
entire alt, tude range. It should be noted that in this format 
distributions occurring in relatively small areas of local time’ 
or latitude w,l| be averaged out in the sort over all local times 
and latitudes above 60°; this is especially true for the M/Z= 2 
ion in which the overwhelming number of nothing observed 
samples masks the contribution of auroral zone outflows 


4.2. Low-Altitude (R E < 3.0) Observations 

Plate I shows low-altitude (R E < 3.0) ILAT versus MLT oc- 
currence probabilities (2.5° by I -hour bins) for rammed (Plate 
la), field-aligned beams (Plate lb), asymmetric (Plate Ic) and 
intermediate distributions (Plate Id). Results are shown for 
H\ He-, and O- and divided into Kp < 3 and Kp > 3. The linear 
scale at the top of these plots and the ones that follow repre- 
sent the occurrence probability as a percentage, the number of 
samples fitting the distribution criteria divided by the total 
number of samples for that ion. The color red represents all 
probabilities greater than the number noted to the right simi- 
larly. the color blue represents all probabilities lower than the 
number noted to the left. Corresponding statistics for M/Z=2 
and O— . without Kp division, are shown in Plate 2. Fewer 
M/Z=2 and 0~ samples at low altitudes preclude statistical cal- 
culations by activity level. Note that, for this altitude range, 
the predawn local time quadrant between 60° and 75° ILAT is 
not well sampled ("white" indicates less than 15 samples were 
obtained in a given bin). Although the discussion applies 
only to the variability of distributions occurring above 60° 
ILAT. we chose 50° ILAT for the minimum value in these spec- 
trograms to provide overlap and ease of comparison with mid- 
latitude statistics reported elsewhere. 

For R f < 3.0 there are relatively few occurrences of H- 
rammed (Plate la) at polar cap or auroral zone latitudes except 
lor a small increase (<20%) in occurrences between 70°-80° 
and 6-12 MLT. All ion species display high occurrences for 
the rammed category below 60° to 65° ILAT. These are plas- 
masphere isotropic distributions and characteristically tend 
toward lower invariant latitudes with increased activity He- 
rammed PADs occur at greater than 50% probability from mid- 
night to dawn between 60° and 80° ILAT. The O- rammed oc- 
curs at >50% probability covering ILAT > 70° from pre 
midnight to about 9 MLT. In addition, the O- rammed occurs at 
>25% probability between 60° and 70° ILAT across the day- 
side. For all ions there is an absence of occurrences at auroral 
zone latitudes for afternoon local times; occurrences also de- 


crease with increasing magnetic activity. Statistics for 
rammed O- follow the O- pattern with the level of occurrence 
probability being similar. WZ=2 follows the He- pattern 
with probabilities lower by -50%. In the case of M/Z= 2 and 
O and to an extent O-, lower probability at plasmasphere 
latitudes as compared to the H- and He- stat, sties will be re- 
flected as high probability for the same latitudes in the low 
counts category. This is a consequence of both the lower rela- 
tive concentration of these ions within the plasma and instru- 
ment operational characteristics (less frequent sampling re- 
suiting in erratic low-count distributions). 

H- low-altitude field-aligned beams (Plate lb) occur with 
greater than 30% probability in an oval pattern that covers 
nearly all local times and appears to be roughly coincident 
with the auroral oval. There is an expansion of this oval to 
lower latitudes for high Kp level and evidence of an increase in 
probability. The He- and O- probabilities also follow the au- 
roral oval pattern except for an absence of beams in the morn- 
ing sector - the area populated by rammed distributions 
Again, there is an expansion to lower latitudes for higher Kp 
and increased probability. O- statistics for the beam distribu- 
tions are not conclusive, although it can be said that O-- 
beams occur with about 10% probability at the same latitudes 
and local times as for O-. Since O- beam statistics increase 
uniformly with activity, we might expect the averaging of 0~ 
beams over all Kp to lower the mean proportions shown in 
Plate 2. The M/Z=2 beams, a category which also contains 
measurements of -1 keV protons, follows the auroral oval 
with -10% probability. 

For all ions the asymmetric PAD (Plate Ic, and also referred 
to as upwelling ion events) has greatest occurrence on the 
prenoon dayside at invariant latitudes above 65° for the light 
ions and above 75° for the heavier ions. In addition, there are 
significant probabilities, increasing with Kp, in the afternoon 
aystde, and, for the light ions, in the region about 20-22 
MLT above 70° ILAT. The overall level of occurrence is 
higher during active conditions by as much as 20%: much of 
this increase being the increased afternoon davside events 
Unlike H- and He-, few O-. M/Z=2. or 0« asymmetric events 
arc observed for the polar cap. Lockwood et al. [1985a] noted 
that when the dayside cleft upwelling feature is identified for 
o , the feature is present for all ions measured, so that the cleft 
region behavior exhibited by O- is representative of that for 
H . He-. M/Z=2, and O- as well. The reverse is not always 
true. When sorting the lighter ions, the programmed pattern 
recognition algorithm identifies additional periods of time as 
fitting the asymmetric criteria (skew coefficient >0.75) 
Therefore more samples and spatial coverage are tagged as 
asymmetric for H- and He- than for O-. The upwelling ion 
events associated with the dayside cleft region [Lockwood et 
ai. 1985a] are clearly seen in the O- statistics but are imbed- 
ded within a larger population of heat flux distributions for H- 
and He + . 

H intermediate distributions (Plate Id) occur at 10-15% 
probability across the dayside between 60° and 70° ILAT and at 
10-15% probability over the polar cap. He- and O- intermedi- 
ate^ are at 10-20% probability over the polar cap and between 
65 to 70° ILAT for all local times. Those over the polar cap 
increase in probability with activity. Statistics for O*- fol- 
lows the O- pattern with a similar spatial pattern but somewhat 
lower probability. Few M/Z= 2 intermediate events are ob- 
served above 55° ILAT. 

Plate 3 shows low-altitude (R E <_ 3.0) ILAT versus MLT oc- 
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eurrencc probabilities for downflowing distributions (Plate 3a) 
and field-aligned conics (Plate 3b). Results are shown only for 
H*. He*, and O* and divided into Kp < 3 and Kp > 3. The 
downflowing distributions (Plate 3a) are difficult to distin- 
guish from the rammed category for the RIMS energy range 
because the downward drift speed of the ions tends to be rela- 
tively small. H* down Hows are discernible at polar cap lati- 
tudes above 80° with an increase in occurrences with activity. 
There are also H* downflows at 70°-75° ILAT from 11 to 13 
MLT. O* downflows appear over the nightside polar cap above 
80° ILAT with no resolvable change with activity increase. 
Irregular regions of O* downflows also appear below 60° ILAT 
and are believed associated with convection near the plasma- 
pause. Occurrences of He*, A//Z=2, and O** downflows are less 
than 5% overall. The H* unidirectional conics (Plate 3b) are 
most frequently observed at invariant latitudes between 70° 
and 80° in the dayside local time sectors, and at invariant lati- 
tudes between 60° and 70° on the nightside. There is an in- 
crease (about 10%) in probability for observations of conics 
with an increase in magnetic activity. The unidirectional 
conic occurrences for the remaining ions are below 10% al- 
though, at least in the case ot He*, the occurrences show hints 
of following the H* pattern. 

4.3. Higher- Altitude (3.0 < R g <^4.67) Observations 

Plate 4 shows high-altitude (3.0 < R E < 4.67) ILAT versus 
MLT occurrence probabilities for rammed (Plate 4a), field- 
aligned beams (Plate 4b), asymmetric (Plate 4c), and interme- 
diate distributions (Plate 4d). Results are for H*, He\ and O* 
ana divided into Kp < 3 and Kp > 3. Corresponding statistics 
for M/Z - 2 and O**, without the Kp division, are shown in Plate 
2 . 

The high occurrences of heavy ion rammed distributions 
(Plate 4a) at low altitudes do not appear in the high-altitude 
statistics. In fact, there is almost a complete absence of 
rammed distributions at auroral zone and polar cap latitudes for 
all the ions measured with the exception of a small area of 
about 10-15% probability for high activity at 70°-80° ILAT 
postmidnight for H + and at 65°-75° ILAT premidnight for O*. 
Again, the probabilities below 60° to 65° ILAT are associated 
with plasmasphere isotropic distributions. 

H + high-altitude beams (Plate 4b) occur not only for auroral 
latitudes but at all latitudes greater than 70° across the full day- 
side sector. On the nightside, beams are observed with compa- 
rable probability only for auroral oval latitudes from midnight 
to dawn, the beams almost disappearing at from dusk to mid- 
night. He* beams occur with greatest probability (20-30%) at 
invariant latitudes between 60° and 80° between 6 to 19 MLT. 
There are somewhat lower, but still significant (15-18%), 
probabilities at 60° to 65° latitude from midnight to dawn. O* 
high-altitude beams occur with greatest probability during 
higher activity. These are at latitudes greater than 70°-75°, 
covering the entire polar cap between for all local times; and 
also from 60° to 65° between midnight and noon. For the 
high-altitude observations, in addition to an invariant latitude 
expansion with increasing magnetic activity, the beams are 
also observed at lower invariant latitudes for heavier masses. 
This expansion is between 5° and 10° for each of the mass 
species. For Af/Z=2, a category which also contains 
measurements of -1 keV protons, there is about 15% 
probability between 75° and 80° ILAT on the dayside and 
between 60° to 70° ILAT on the nightside. The probability for 



Occurrence Frequency 



field-aligned beams, (c) asymmetric, and (d) intermediate distributions. Results are for H + , He + , and 0 + and 
divided into Kp < 3 and Kp > 3. 
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°** unidirectional beams is less than 10% with the exception 
of a small area of about 15% probability at 70°-80° ILAT and 
11-13 MLT. 

The asymmetric distributions (Plate 4c) make up less than 
5% of the PAD samples at high altitude with the exception of a 
small area ot about 10% probability for H* at 65°-75° ILAT on 
the nightside for high magnetic activity. Similar to the ram 
and asymmetric results, the intermediate category (Plate 4d) is 
a much smaller proportion of the high-altitude sample popula- 
tion and, instead of appea'ing spread over the entire polar cap 
as tor the low-altitude observations, they are limited to the 
nightside auroral zone between 20 MLT and 2-4 MLT, The dis- 
tribution of observations is similar for all ions with lighter 
ions appearing at slightly lower latitudes and movement to- 
ward higher latitudes with higher magnetic activity, 

4.4. Absence of Cold Plasma Observations 

Plate 5 shows low-altitude (fi £ < 3.0) ILAT versus MLT oc- 
currence probabilities for nothing observed (Plate 5a) and low 
count distributions (Plate 5b) and high-altitude (3.0 < R £ < 

4 67) probabilities for nothing observed (Plate 5c) and low 
count distributions (Plate 5d). Results are shown for H*. He*, 
and O* and divided into Kp < 3 and Kp > 3. Statistics for 
,V//Z=2 and O**, not sorted for activity, are shown in Plate 2. 

H* low-altitude nothing observed events (Plate 5a) occur for 
latitudes above 70° ILAT across the nightside and above 80° 
ILAT across the dayside for low Kp and with lower frequency 
above 70° ILAT at all local times for high Kp. For high alti- 
tude the H* nothing observed events (Plate 5b) occur over the 
polar cap and, for low Kp, at auroral zone latitudes near mid- 
night. There is also significant probability below 60° ILAT at 
dusk, indicating a depletion of thermal ions in the plasma 
trough area there. He* low-altitude nothing observed events 
occur at >30% probability for all latitudes above 65° across 
the nightside and above 70° ILAT across the dayside. There is 
an expansion to slightly lower latitudes for higher Kp and 
slightly overall decreased probability. For high altitude the 
He* nothing observed events occur at >50% probability for all 
latitudes above 65°-70° ILAT at all local times. High-altitude 
O* nothing observed events occur at up to 75% probability 
over the entire polar cap and auroral zone area for low Kp\ for 
high Kp the O* nothing observed events drop dramatically in 
frequency and follow the auroral oval pattern. O** generally 
follows the pattern of O*. although without division with Kp it 
is not known whether a decrease in occurrence probability oc- 
curs with activity. M/Z= 2 follows that of He*. 

Low-count statistics are not significant at low altitude 
(Plate 5c), the observations that do occur appear at outer plas- 
masphere/plasma trough latitudes. For high altitude the low 
count category (Plate 5d) becomes important for the H* ion 
near midnight at auroral zone and polar cap latitudes. 

5. Discussion 

The low-energy (<50 eV) H*. He*, O*, M/Z= 2, and O** RIMS 
spin angle distributions for latitudes greater than 60° have 
been classified into the following pitch angle distribution cat- 
egories. unidirectional field-aligned beams, downllows, con- 
ics. asymmetric, rammed, intermediate, low counts, and noth- 
ing observed. 


5.1. Ionospheric Source Regions 

The statistical survey shows that low-enerey outflow from 
the high-latitude ionosphere occurs for all the measured ions 
with a spatial and statistical distribution dependent on mae- 
net’c activity. For convenience, we discuss the outflows in 
terms ol three regions: the dayside cleft, auroral zone, and po- 
lar cap. 

During quiet times, cleft region upwelline ion events (the 
asymmetric distributions in Plate 1) are observed almost ex- 
clusively in the prenoon sector. The area of most probable 
observation becomes smaller and is bounded at higher invari- 
ant latitudes for the heavier mass species. These asymmetric, 
field-aligned spin distributions are identified as ions carrying 
an upward heat flux, with both parallel and perpendicular ve- 
locity components, and so "upwell" in response to heating be- 
low the satellite [Moore et al. , 1985; Biddle et al. , 1985; 
Lockwood et al., 1985a]. Moore et al. [1985] demonstrated 
that the upweilings exhibit conic structure in the energetic tail 
of the distribution and suggested that the distinction between 
the unidirectional beam, conic, and upwelling categories may 
be a matter of the degree by which transverse acceleration af- 
teas the low-energy core of the distribution. The probability 
of observing upwelling events in the morning sector is rela- 
tively steady regardless of activity level which contrasts 
sharply with the large increase with activity seen with the au- 
roral CT outflowing beams (Plate lb). This is consistent with 
Lockwood et al. [1985a] in which these prenoon 0 + outflows 
from the polar cleft ionosphere were first identified as a major 
source of heavy ions for the polar cap magnetosphere. There 
is, in addition, a distinct spatial change for cleft outflow with 
increased activity - at higher Kp levels upwelling events ap- 
pear on the afternoon dayside and the area of peak probability 
spreads to lower latitudes. In a study of the source characteris- 
tics of these events, Moore et al. [1986] found typical ener- 
gies to be of the order of 10 eV and found that the distinctive 
transverse ion heating exceeds gravitational escape require- 
ments. Evidence from this and the other studies indicates 
then, that the dayside prenoon cleft is a consistent source of 
low-energy ions of velocity sufficient to escape to outer mag- 
netosphere regions regardless of species (including the minor 
ions M/Z=2 and 0 + +) or activity level and that the dayside af- 
ternoon auroral zone becomes an additional source during in- 
creasing activity. 

While it is clear from the statistics that this asymmetric 
data signature, an indicator of transverse heating in the iono- 
spheric outflows, is specific to auroral latitudes, it cannot be 
unambiguously stated that they are unique to the dayside 
auroral zone. Figure 6b shows the peak observance of these 
events to be between 1.5 and 2.0 R E geocentric altitude. 
Unfortunately, due to a combination of instrument operation 
and orbital constraints, samples taken below 1.75 R E geocen- 
tric cover only the 5 to 12 MLT range. Samples between 1.75 
and 2.0 R E geocentric add auroral latitudes from 19 to 24 MLT. 

If, as has been suggested, the upwelling data signature evolves 
by magnetic folding to a field-aligned conic distribution, and 
if the asymmetric data signature does exist for the entire auro- 
ral zone, we would expect to find ample evidence of conics at 
nightside auroral latitudes, which is the case as can be seen in 
Plate 3b. 

Low-energy ion outflows associated with the auroral oval it- 
self are mainly of the unidirectional beam type (Plate lb). For 
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H . unidirectional beams are observed in the auroral zone tor 
most ail local times, increasing by about 10% with activity. 
The beam distributions arc observed at progressively lower 
invariant latitudes for heavier mass species. These 
observations are in agreement with the earlier studies of 
Gluclmem et ai [1978] and Gome y et ai [1981 ] in that the 
occurrence ol upward flowing beams appears to be associated 
ttdth the statistical auroral oval. The occurrence frequencies 
peak more equatorward on the mghtside than for the dayside. 
and peak at lower invariant latitudes during disturbed times. 
The observations are also consistent with Yau et oi [1984], 
who also found that the oval of upflowing beams is energy de- 
pendent (the iow-energy upflowing oval occurring poleward of 
the energetic ion upflowing oval). In addition, for heavier 
ions the beams display a distinct asymmetry in distribution 
speed between dawn and dusk local time sectors. For He* and 
O the unidirectional beams are observed primarily for the dusk 
sector and increase significantly in probability with increas- 
ing activity (as much as a 30% increase for 0 + ). The "nothing 
observed” category complements this duskside behavior by 
having probability greater than 50% in this same region for 
low activity and decreasing sharply for high activity (Plate 
.''a). In the dawn sector the heavier ion distributions, espe- 
cially O**. tend to be the lower velocity rammed type (Plate 
la). Since the unidirectional beam distributions generally 
have adequate escape velocity (Figure 4), the auroral zone then 
appears to be a steady source of escape velocity H* for all local 
times, a steady source of escape velocity He* ions on the 
duskside, and a source of escape velocity heavy ions on the 
duskside but primarily during increased activity. The fact that 
the He + and O* low-altitude dawn-sector rammed and 
intermediate distributions (Plates la and Id) fill in the missing 
beam distributions suggests that there is an amount of low- 
energy heavy ions unable to escape to high altitudes unless 
acted upon by energization processes which are correlated with 
increased activity. The M/Z=2 ion beam statistics contain 
measurements of ~1 keV protons and so it is unclear whether 
escape velocity outflows other than the asymmetric type exist 
lor this mass type. 

l or latitudes poleward of the auroral zone the pitch angle 
distribution statistics vary with mass, activity level, and local 
time. The majority of the heavy ion observations in which 
significant counts are detected (>50 or 60%) are of the rammed 
and intermediate type (Plates la and Id), which still indicate 
I ield-al igned upflow, but with velocities comparable to the 
spacecraft and shown in section 3 not to be sufficient for es- 
cape to higher altitudes without additional energization. There 
is also a dawn-dusk asymmetry reflected in the ion distribution 
’speed" in which higher speed distributions appear with great- 
est frequency at dusk and lower speed distributions have great- 
est frequency at dawn. There appear to be small changes in the 
probability levels with magnetic activity, a slight increase in 
rammed distributions, and a small decrease in intermediate ob- 
servations; however, it is important to note that there is al- 
ways low-energy O* at low altitudes over the entire high- 
latitude region regardless of activity (Figure 6b). On the other 
hand. H* has a significant number of heat flux (asymmetric), 
unidirectional beams, and intermediate events for latitudes 
poleward of the auroral zone, the probabilities for these cate- 
gories increasing with activity. There are relatively few 
rammed distributions H + , although the "nothing observed" 
category for H + has >50% probability for low Kp. The He* ob- 
servations have characteristics of both H* and O*, the "noth- 


ing observed /heat flux statistics of H* and the 
rammed/intermediate statistics of O*. Overall, the region 
poleward of the auroral zone is again a consistent source of 
low-energy tons, although we would expect only the lighter 
mass particles to have sufficient velocity to escape to higher 
altitudes. 

5.2. Evolution of Field-Aligned Distributions 
During Outflow 

By comparing distribution types at low altitudes with those 
at high altitudes we can infer something about how low-energy 
ion distributions evolve after ionosphere injection. The 
rammed distributions, which make up a significant portion of 
dawnside heavy ion observations at low altitude, essentially 
disappear at high altitude (Plate 4a). These either reenter the 
ionosphere or are trapped by magnetic mirror effects. If these 
low drift velocity rammed distributions instead experienced 
energization, such that the ramming effect no longer domi- 
nated, then we would expect a significant portion of "faster 
speed" distributions such as intermediate or field-aligned 
beams in the dawnside polar cap, which is not the case. For 
similar reasons the intermediate distributions seen at low alti- 
tudes are also believed to return to the ionosphere or remain 
trapped. At low altitude, only the asymmetric and field- 
aligned conic distributions, which reflect a significant degree 
of transverse acceleration, and the field-aligned beam distribu- 
tions appear to account for the statistical patterns observed at 
high altitudes. The observations support two mechanisms, 
which may well be related, for low-energy outflow into the 
outer magnetosphere (the cleft ion fountain, so named by 
Lockwood et ai. [1985b], and its counterpan in other sectors 
of the auroral oval) what the present authors call the auroral 
ion fountain. 

The cleft ion fountain, as explained by Lockwood et al. 
[1985b], consists of ionospheric plasma of 1-20 eV energy 
streaming upward into the magnetosphere where high-latitude 
convection electric fields can modify the ion path resulting in 
a poleward dispersion. This effect is a well-documented feature 
[cf. Waite et al. . 1986; Moore et al., 1985; Tsunoda et ai, 
1989] and is clearly evident by the beam distributions at high 
altitude over the polar cap (Plate 4b). It is not certain whether 
all the high-altitude O* polar cap observations originate from 
the cleft ion fountain mechanism; however, the very low 
probabilities of observing escape velocity O* outflow at low 
altitudes in the polar cap would indicate that this is so. The 
outflows in Plate 4b also demonstrate the mass spectrometer 
effect associated with the cleft fountain whereby the lighter 
ions escape to the magnetosphere tail regions, but the heavier 
low-energy ions, such as O*. are more influenced by the con- 
vection electric field and are swept toward the nightside auroral 
region and plasma sheet. For weaker convection, heavier ions 
outflowing on the dayside would be swept antisunward over the 
polar cap and reenter the nightside polar ionosphere. This is 
supported by the Plate 3 O* downflows. 

The aurora] ion fountain concept is reflected in the field- 
aligned beam and conic results. For the RIMS instrument the 
unidirectional field-aligned beams correspond to pitch angle 
distributions that have not undergone magnetic mirroring and 
can be represented by a drifting Maxwellian distribution with 
supersonic drift speed [ Moore et ai, 1985]. Traveling upward, 
these ions flow along high-latitude field lines that are open; or 
closed in a highly elongated fashion such that bouncing would 
occur over long time periods and be disrupted by nonadiabatic 
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processes. Low-altitude auroral latitudes are populated by 
beam distributions near dusk, "slower" intermediate or rammed 
distributions at dawn, and conics on the nightside (Plates 1 
and 3). High-altitude, late-aftemoon beams are almost absent; 
instead the region is represented by the intermediate and 
"nothing observed" categories (Plates 4 and 5). It is hypothe- 
sized that low-altitude field-aligned beams flow outward along 
auroral field lines, and for late afternoon, experience addi- 
tional acceleration, such that the ion energy distribution ex- 
ceeds the instrument detection range (>50-60 eV). Higher- 
energy beams, particularly those >1 keV, have been shown to 
have strong local time dependence, being concentrated near 
dusk at auroral latitudes [Gomey et ai , 1981; Yau et ai , 1984; 
Kondo et ai , 1990], 

The sequence of events for both the cleft ion fountain and 
the auroral ion fountain is consistent with the theory that a 
continuum of behavior exists such that the upwelling distribu- 
tions evolve with altitude into conics and then high-altitude, 
t ield-al igned flows (cf. Figure 6). Conics correspond to up- 
ward field-aligned flows that have some additional perpendicu- 
lar velocity component. Conic distributions typical of the 
midaltitude nightside have been studied by Moore et ai [1985] 
and found to originate between 1000 and 2000 km altitude 
which is consistent with the altitude distribution shown in 
Figure 6. The study of Gomey et al. [1981], based on electro- 
static analyzer data at higher energies (0.09 to 3.9 keV), found 
that ion conics occur in the geocentric distance range of 1.4 up 
to the S3- 3 apogee at 2.3 /? f , or about 2500- to 8000-km alti- 
tude and that ion beams occur primarily above about 2.2 R E or 
7200-km altitude, with the frequency increasing with altitude. 
During quiet times the S3-3 conics were associated with the 
daytime polar cusp region, but were observed uniformly in lo- 
cal time during disturbed times. The older studies did not dis- 
tinguish between events peaked near 90° pitch angle (our 
asymmetric category) and those peaked between 120° and 170° 
pitch angle (our conic category). However, they are consis- 
tent with our observations in that our dayside asymmetric 
events exist for all activity levels, while those observations 
associated with our more narrowly defined conic category are 
found to increase with activity and are distributed in local 
time. The S3-3 results are also consistent in that quiet time 
ion beams have maximum frequency in the premidnight sector 
and are mainly confined to the dusk sector during disturbed 
times. 

5,3. Considerations on the Absence of Cold 
Plasma Observations 

The "low-count" and "nothing observed" categories are rep- 
resentative of two possible conditions: regions in which the 
ions have energy so low that they are repelled by positive 
spacecraft potentials or of regions where there are no ions 
within the detection range of the instrument, the low-count 
samples being those that occur near the limits of the instru- 
ment response (-0.1 ions enr 3 with no potential barrier). 
Analysis of RIMS data by Olsen [1989] established a density- 
potential relationship such that the satellite accumulates posi- 
tive charge on entering density regions below 1000 cm 3 , ris- 
ing slowly to about +1 V at 100 cm \ and about +5 V at 10 
cm 3 . By preventing measurements of the coldest plasma 
components, this potential can mask out an isotropic 
background plasma or, when the potential is sufficiently high, 
can give the appearance that no cold plasma exists in a region. 


This effect is known to be a concern in the plasmapause and 
trough region, where densities are of the order of 7 to 70 cm 3 
and the spacecraft reaches potentials of +2 to +5 V [Olsen, 
1989]. At altitudes above the ionosphere in the polar regions, 
densities again can be sufficiently low so that observation of 
ions of classical polar wind speed is possible only during 
negative aperture bias instrument modes to offset the positive 
potential [Nagai et ai , 1984]. We cannot always assume, in 
low-density regions, that cold plasma does not exist just 
because it has not been measured. 

The RIMS H* nothing observed events at low altitude may 
be representative of an inability to observe the lowest-energy 
polar wind outflows in low-density regions. As detailed in the 
above discussion on the "auroral ion fountain," at high alti- 
tude near midnight, the H* nothing observed events are 
thought to be measurements in which the low-altitude outflows 
have been accelerated beyond the range of the instrument. The 
He* and 0 + events in the dusk auroral zone are probably a mix- 
ture of the two effects. 

6. Summary 

We have presented a statistical study of the occurrence fre- 
quency distribution of low-energy ion pitch angie distribu- 
tions in the near-Earth magnetosphere using measurements 
from the retarding ion mass spectrometer aboard Dynamics 
Explorer 1. A relationship is shown between pitch angle dis- 
tributions of different masses (H + , He*, O*, M/Z-2, and O**) 
and velocity characteristics. The added spatial dimensions in 
this study relative to earlier work provide a more detailed pic- 
ture of ion transport events. Occurrence frequencies are deter- 
mined for all measured ions as a function of altitude, invariant 
latitude, magnetic local time, and magnetic activity (Kp). 
From this study, several characteristics are apparent. 

6.1. Overall Occurrence With Activity 

The probability of observing H* and M/Z- 2 distributions 
tor 1LAT >60° was found to be independent of magnetic activ- 
ity. He* was found to be weakly dependent, while O* and O** 
were strongly dependent on magnetic activity. The large in- 
crease in O* and O** observations with increasing activity was 
found to be composed primarily of escape velocity unidirec- 
tional Field-aligned and upwelling distribution types. 

Upwelling (transversely heated, upflowing) ions occur pre- 
dominantly, but perhaps not exclusively, in the pre-noon sec- 
tor at 70°-80° ILAT. At higher Kp levels the events are ob- 
served with increased probability on the afternoon dayside and 
at slightly lower latitudes for all ion species. 

The occurrence of unidirectional beams, for all ion species 
including M/Z- 2, is consistent with the location and geomag- 
netic activity dependent expansion and contraction of the au- 
roral oval. The data suggest that beam distributions observed 
at higher altitudes are predominately associated with the up- 
welling distributions observed at lower altitudes and that the 
paths for ions injected at low altitude in the dayside auroral 
oval include some amount of convection antisunward over the 
polar cap, as suggested in the Lockwood et al. [1985a, b] cleft 
ion fountain sequence. 

Peak conic and upwelling ion observations occur predomi- 
nantly at geocentric distances of 1.75 and 2.0 R E , respec- 
tively, while field-aligned beams occur most frequently above 
about 2.0 to 2.25 R E , especially for the heavier ions. Since 
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angular distributions tend to become increasingly "beamlike" 
at higher altitudes because of magnetic focusing it is thought 
that the upwelling events may well evolve into the observed 
unidirectional conic and beam distributions as the ions travel 
outward along the magnetic field direction. 

6.2. Ionospheric Source Regions 

The dayside prenoon cleft is a consistent source of low- 
energy ions of escape velocity regardless of species or 
activity level and the dayside afternoon cleft, or auroral zone, 
region adds to this source for increased activity. 

The auroral oval as a whole appears to be a steady source of 
escape velocity H + ions at all local times, a steady source of 
escape velocity He + ions for the dusk sector, and a source of es- 
cape velocity heavy ions for dusk but primarily for increased 
activity. 

The region poleward of the auroral zone is a consistent 
source of low-energy ions, although only the lighter mass par- 
ticles appear to have sufficient velocity to escape to higher al- 
titudes. 

6.3. Evolution of Field-Aligned Distributions 
During Outflow 

The observations support two concepts for outflow: (1) The 
cleft ion fountain consists of ionospheric plasma of 1-20 eV 
energy streaming upward into the magnetosphere where high- 
latitude convection electric fields modify the ion paths result- 
ing in a poleward dispersion of the ion trajectories. (2) The 
auroral ion fountain consists of field-aligned beams which 
stream outward along auroral latitude field lines, and in addi- 
tion, for late afternoon local times, experience additional ac- 
celeration, such that ion energy distributions tend to exceed 
the detection range of the instrument (>50-60 eV). 

6.4. Absence of Cold Plasma Observations 

The low-count and nothing observed’ trends with mag- 
netic activity appear to complement that of the unidirectional 
field-aligned beam distributions and are representative of two 
possible conditions: regions in which the ions have energies 
relative to the spacecraft so low that they are repelled by posi- 
tive spacecraft potentials or regions with no ions having en- 
ergies within the detection range of the instrument. We can- 
not always assume, in low-density regions, that cold plasma 
does not exist just because it has not been measured. 
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The characteristics ot the polar ion outflows as observed in the topside polar ionosphere by the Dynamics Explorer 
Retarding Ion Mass Spectrometer are reported in this study. The study is restricted to altitudes between 1000 and 
4000 km in order to locus on the phenomenon ol the "classical polar wind.” Except insofar as they convect into the 
polar cap at these altitudes, auroral and cusp/cleft outflows are specifically excluded from this study. Using a 
method based on the "relative wind" ot ions as seen from the moving spacecraft, averages and variances of the 
magnetic field-aligned ion flux and velocity, and the species densities, are derived and binned for examination of 
their altitudinal, seasonal, and magnetic activity dependencies. The data set used encompasses the period from 
autumn 1981 through 1983. the altitude range trom 1000 to 4000 km. and invariant latitudes greater than 70°. 
Observations ot . He“. and O* with energies down to spacecraft potential are reported. It is found that the 
polar outflows are not well correlated with common indicators of solar or geomagnetic activity. The average 
parameters ot the outflows exhibit a winter enhancement of He'*' flux, a winter decrease in H* flux, and a 
summer enhancement ot the asymptotic outflow speeds of all species. The solar F , 0 7 (and hence EUV) flux is 
correlated with the acceleration profile of the flows, with slower flow velocities at 2500 km altitude for active solar 
conditions. Evidence is found that the H~ flows are. at times, subsonic in the altitude range studied, in contrast to 
cold polar wind models. 


Introduction 

The upward flow of ionospheric ions from Earth’ s polar region 
has been studied extensively over the past two decades with the 
bulk of the work being theoretical. The earliest suggestions regard- 
ing this phenomenon were based on the then new idea of open 
magnetic field lines connecting the polar ionosphere to the distant 
magnetotail and the solar wind magnetic field which would provide 
a means of escape for ions energetic enough to overcome Earth’ s 
gravitational pull. Dungey 1 1 96 1 ] noted that open magnetic field 
lines should upset the plasma equilibrium along those field lines. 
Dessler and Michel 1 1966] and Nishida ( 1 966) pointed out that at 
sufficiently high temperatures H + would “evaporate" from the 
ionosphere and escape into the tail. Numerical studies involving 
solutions to the steady state hydrodynamic equations showed that 
when the dominant species in the ionosphere is O* these heavier 
ions and electrons play a role in the ion outflow process by creating 
an ambipolar electric field through which the lighter H * and He + 
arc accelerated and Hoher, 1 968, 1969 a,b\. Axford 1 1968] 
invoked this “polar wind’’ flow to provide a nonthermal escape 
mechanism for He from the terrestrial atmosphere. Following these 
initial modeling efforts, discussions of the outflow under collision- 
's conditions spawned several kinetic models | Dessler and 
Cloutier, 1969: Letnaire and Scherer, 1970. 1971. 1972] which 
were later shown to be in agreement with the hydrodynamic formu- 
lations \Banks et al., 1971]. 

Increasingly, more detailed models have been developed to 
s ludy the dynamics and energetics of the outflow of H* , He + . 
and CT from the polar region \Marubashi, 1970; Hoher et al.. 
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1971; Lemaire , 1972; Banks, 1973; Bailey and Moffett, 1974; 

Banks et al., 1974, 1976; Strobe I and Weber, 1972; Raitt et al., 

1975, 1977, 1978^,6; Schunk et al., 1978; Schunk and Watkins, 

1 979, 1981, 1 982: Ottley and Schunk . 1 980: Gombosi et a!., 1 985; 

Barakat et al., 1987; Schunk and Sojka, 1989]. All but the most 
recent of these have been reviewed by Raitt and Schunk [1983]. 

The basic predictions for ion outflow have remained roughly the 
same over two decades at — 1 0 8 and — 1 0 7 ions cm~ : s" 1 (at 1 000 km 
altitude) for H + and He + , respectively. 

Raitt and Schunk |1983| have pointed out the lack of experi- 
mental observations of the polar cap outflows. To date, the results 
include initial observations of the nondiffusive density profiles 
and upward H~ flow by Hoffman |1970| from Explorer 31 and 
measurements of the outflow of H* and He' from ISIS 2 by 
Hoffman et al. |1974| and Hoffman and Dodson 1 1980|. This 
latter study included data from a 4-week period centered on the 
winter and summer solstices and the spring and autumn equinoxes 
of 1971 and 1972 with approximately 50 orbits of data for each 
season at an altitude of 1400 km. The basic results include H* 
fluxes of the order of 10* cm" : s" 1 and He* fluxes of the order of 
10 cnT : s _l (at 1000 km altitude) with strong seasonal 
dependence. 

Case studies of polar outflow have recently been carried out 
using data from the Dynamics Explorer (DE) mission | Gttrgiolo 
and Burch. 1982; Sojka et al.. 1983: Hagai et al.. 1984]. These 
studies concentrated on the ion velocity distributions of polar 
wind outflows demonstrating the supersonic behavior of the H 
ion outflow at high altitudes. Biddle et al. |1985| also used DE 
data to determine the ion heat flux in polar wind outflows and 
found at least one case to be in agreement with theoretical predic- 
tions that the sense of the heat flux should reverse as the plasma 
flow becomes supersonic. In addition, results of statistical studies 
of upflowing ions with energies greater than 10 eV have been 
presented by Yau et al. |l985«.fr). 
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The following results arc based on observations of polar cap ion 
outflows using the retarding ion mass spectrometer (RIMS) on DE 
1 \C lutppell et aL. 19811. The purpose is to show average values 
of ion flux, velocity, and density as well as altitudinal, seasonal, 
and magnetic activity-related variations in the polar outflows for 
energies down to ionospheric values (<1 eV). The addition of 
these results to those already published will provide a more 
complete picture of the phenomena of ion outflow from the polar 
region. 


Data and Techniques 

The data for this study were taken during the period from 
October 1981 through April 1983 and cover the altitude range 
from 1000 to 4000 km. All observations are from the RIMS radial 
detector which samples ions in the satellite spin plane. The spin 
plane contains the satellite velocity vector and is parallel to the 
orbit plane. The ion velocity and flux measurements are obtained 
by integrating the observed flux over the spin cycle of the satellite 
for a zero retarding potential analyzer (RPA) voltage. A full 
description of the technique with error analysis can be found in 
the work by Chandler and Chappell [1986]. 

It should be noted that using the zero RPA setting implies that 
any superthcrmal population which may be present in the high- 
latitude ionosphere will also be included in this analysis. This is 
unavoidable since these data are taken from the period when the 
RPA was not functioning. An attempt has been made to minimize 
the possibility of “contamination*' by superthermal populations of 
either auroral or cusp origin by restricting the observations to 
invariant latitudes greater than 70° and altitudes below 4000 km. 

The velocity measurements were made under the assumption of 
zero convective plasma motion in the satellite spin plane. Error 
analysis shows that this assumption results in errors of less than 
50# in the resulting velocities for convective electric fields, at 
ionospheric altitudes, of less than 100 mV m To minimize the 
effects of convective electric fields the data have been restricted to 

invariant latitudes greater than 70°. 

For the average flux and velocity values the statistical standard 
deviation is typically of the order of the average value of the quan- 
tity. Thus all error bars shown in this paper are the statistical 
values for one standard deviation. Chandler and Chappell [ 1986] 
concluded that the error in determining the ion velocity was of the 
order of 100 m s' 1 . In the present study the data have been aver- 
aged into 2°-wide bins in spin-phase angle. This ± 1° uncertainty 
in the angles yields errors of the order of 20# following the error 
analysis of Chandler and Chappell [1986]. However, it also 
yields a lower limit on the measurable velocity of -150 m s' 1 
(under the assumption that the satellite velocity is 8 km s 1 and 
perpendicular to the magnetic field). Given the high outflow 
velocities for H* and He* the 20# error is the more useful 
number. For O*, however, this lower limit results in some 
ambiguity for the lowest measured values. 

The calculation of ion densities in this study assumes zero 
spacecraft potential. Because of the large spacecraft velocity the 
O* ram energy was >4.5 eV in all cases. Thus the O* Mach 
number, relative to the spacecraft, was >2 (assuming an O* 
thermal energy less than 1 eV) and only small errors (<20#) in 
the 0 + density calculations are expected due to spacecraft poten- 
tial. Likewise, the He* Mach number was likely to be >1 result- 
ing in errors less than a factor of 2 for potentials less than 2 V . On 
the other hand, the H * Mach number could range from I down to 
0.5 for a I eV thermal energy resulting in density estimates too 
low by up to a factor of 4 for potentials up to 2 V. However, the 


average O* density in this data set is >500 cm 3 at all altitudes 
which, based on the results of Olsen [1989], gives an estimated 
spacecraft potential of less than 1 V. For a I V potential the max- 
imum errors for He* and H * densities become factors of 1 .3 and 
2.0, respectively. These estimates of the reduction of flux due to 
positive spacecraft potential are derived from simulations of the 
RIMS instrument using the equations of Parker and Whipple 
[1970]. 

Approximately 2000 l-min-duration observations were 
analyzed for each ion. After restricting the data to high latitudes 
the total number of observations is 914 for each ion species. The 
distribution of data samples with respect to season was 61 samples 
in spring, 205 samples in summer, 521 samples in autumn, and 
127 samples in winter. 

Observational Results 

Average Flux 

Variability with geophysical parameters. It is expected that the 
outflow of ions from the polar regions will vary in its magnitude 
with variations in geophysical parameters. The RIMS polar wind 
data have been studied in an attempt to show the magnitude of 
such variations in velocity and flux. The polar wind outflow 
should vary with season as a result of changes in the solar EUV 
input as well as from known seasonal variations in the neutral 
atmosphere composition (e.g., winter enhancements in [H] and 
[He]). Table 1 shows the average flux of H* and He* for four 
seasons defined by contiguous periods centered on the equinoxes 
and the solstices. The H * flux exhibited a factor of 3 variation over 
the seasons, being largest (1.5 x 10 8 cm" 2 s 1 ) in spring and 
smallest (5 X I0 7 cm -2 s' 1 ) in the winter. In contrast, the He* flux 
was observed to reach its maximum value in the winter (2X 10 ^ 
cm -2 s' 1 ) and its minimum value in the summer ( I x 10 6 cm - sVjKg 
In addition, the He* winter value was larger than that observedTiL 
spring by a factor of 2. 


TABLE I . Seasonal Variation of the Polar Wind Flux 


Season 

Flux, cm* s 1 



He + 

Winter 

6x I0 7 

2 x 10 7 

Spring 

l.5x 10* 

1 x 10 7 

Summer 

5 x 10 7 

1 x I0 6 

Autumn 

8x I0 7 

6x I0 6 


Solar cycle variations. The RIMS data in this study come from 
a period of declining solar activity. However, because of the rela- 
tively short time period spanned by the observations, the variation 
in the solar 10.7 cm flux was not large. There are few observa- 
tions from periods of F ]0 .7 > 200 (1CT 22 W trf 2 Hz' 1 ) and fewer 
still for F10.7 < 100 (10 -22 W nrf 2 Hz' 1 ). Also, for the winter por- 
tion of the year there was insufficient spread in the F | 0 7 values to 
make meaningful correlations with the magnitude of the flux. For 
the summer season the spread was large enough to provide results, 
however, there is no indication of a significant trend in the data 
with respect to F ]0 j- 

Variation with magnetic activity'. To determine the effects of 
geomagnetic activity on the ion flux in the polar region, two 
indices were chosen, Kp and Ap. In the correlations of flux with 
the Kp index it was found that using the index from the 3-hour 
period previous to the time of a given observation resulted in ave *| 
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ages with smaller standard deviations. While there exist some 
variations in the observed outflow with the Kp index, the vari- 
ation is only a factor of 3 over the entire range of values observed 
(0 to 7-). Similarly, the variations in H + with the daily-averaged 
Ap also show about a factor of 3 variation. Neither of these vari- 
atiions is considered significant when compared to the standard 
deviations about the means. 

Average Velocity 

The average H* upward velocity, from all data at invariant lati- 
tudes greater than 70°, is shown in Figure 1 . At altitudes below 
-2000 km the average speed is near constant at 3-4 km s“‘ . In the 
region between 2000 and 3000 km the speed increases to — 10 km 
s _l and remains nearly constant, albeit with an increased standard 
deviation, above 3000 km. In order to illustrate the effects of dif- 
ferent geophysical parameters on the ion speed, the data was 
binned according to the day of the year and by the value of the 
10.7-cm solar flux. Figure 2 shows the H * velocity averaged over 
day-of-year for values of the 10.7 cm flux <140 (1CT 22 W m~ 2 
Hz -1 ) and >190 (lO -22 W m -2 Hz" 1 ). While this binning excludes 



1 000 2000 3000 4000 

Altitude (km) 

Fig. 1. The average velocity of as a function of altitude for all data. 



1000 2000 3000 4000 

Altitude (km) 

Fig. 2. The average velocity of H * as a function of altitude for F )0 7 < 140 
(lO 22 W m 2 Hz ') and F ti)7 > 190 (10"" W nr 2 Hz 


a significant number of observations which fall in the range from 
1 40 ( 1 0~ 22 W m " 2 Hz ') to 190 (1CT 22 W nT 2 Hz“‘), the separation 
between the bins accentuates the differences. The difference 
between these two profiles appears in the midaltitude region 
where for low F m 7 the speeds in this region are —509c higher than 
for high F|q 7 flux. In addition, in the high F 10.7 flux case the 
acceleration occurs above 2000 km whereas for the low F ]01 flux 
case the increase in speed begins below 2000 km. While this 
effect is small and is within the standard deviations, the fact that it 
is seen in all the averages between 1500 and 4000 km suggests 
that it is indeed a real effect. 

For averages over all F !0 7 values but at different seasons of the 
year the profiles again exhibit significant differences. In Figure 3 
the speeds are shown for the summer season (May 6-August 4 in 
the northern hemisphere and November 4-February 3 in the 
southern hemisphere) and the winter season (November 
4-February 3 in the northern hemisphere and May 6-August 4 in 
the southern hemisphere). There are no significant differences in 
the average speeds in the altitude range below 2000 km. How- 
ever, at altitudes above 2500 km the average speed for summer is 
—25-70% higher than the winter average. Again, as in the case of 
the results in Figure 2, the deviation about the mean is large in the 
high-altitude region, suggesting additional influences controlling 
the maximum speed. 



1 000 2000 3000 4000 

Altitude (km) 

Fig. 3. The average velocity of H * as a function of altitude for the summer 
season and the winter season. 


The two-season binning is not entirely satisfactory for showing 
the real seasonal variation. Unfortunately, there are insufficient 
data to allow meaningful results as a function of altitude to be 
obtained for all seasons. There are sufficient data to produce good 
statistics for the autumnal equinox case. As might be expected the 
maximum average speeds fall between those for summer and 
winter. 

The outflow speeds for He + are about a factor of 3-4 less than 
the values on average, although the deviations about the mean 
values are similar. While the altitude profiles arc similar to those 
of H " . there exist some noticeable differences. Figure 4 shows the 
He + speeds for high and low F, 0 7 flux. Due to the larger variations 
about the mean values these results arc not as clean as those for 
H + . The results arc. however, similar to those for H ^ with regard 
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Fig. 4. The average velocity of He* as a function of altitude for F io 7 < 
(10-“ W m"' Hz ’) and F, 07 > 190 (I0“ 22 W m“ 2 Hz' 1 ). 


140 


to the F I0 7 variation. As for the seasonal variation (Figure 5) the 
difference between summer and winter is less than 209c (with the 
summer velocities lower than the winter velocities) which is much 
less than the standard deviations. 

In the majority of the observations the 0 + outflow speed was 
less than 2 km s _l so that experimental uncertainties make a 
comparatively larger contribution to the standard deviations. Thus 
the statistical significance is less than for or He + . 

Velocity' correlations . Over the range of altitudes covered by 
this study, 0 + was the dominant ion. Therefore it is expected that 
both H + and He + would behave as minor ions. As the previous 
figures have shown there are similarities in the speed-vs-altitude 
profiles for the two minor species. To investigate the degree of 
correlation between the velocities for the two species, the H + 
speeds for altitudes between 2000 and 4000 km have been sorted 
into bins with 2 km s" 1 widths and averages have been computed 
for the corresponding He* speeds. Figure 6 shows the results 
along with a linear fit to the data. The resulting parameters arc 
given in the figure. 


7 8h 
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Fig. 5. The average velocity of He * as a function of altitude tor the summer 
season and the winter season. 



Fig. 6. The average velocity of He* as a function of the H^ velocity. 


The results in Figure 7 show H + velocity vs. O* velocity for 
altitudes greater than 2000 km. While there is a systematic 
increase in the average H + speed with the change in O* velocity 
from -2 km s _1 (i.e., downward) to + 2 km s‘ l (upward) the effect ^ 
is not striking. However, the range of observed H* speeds, as 
evidenced by the bars corresponding to the standard deviation 
about the mean value, also increases as the O* speed increases 
from -2 to 2 km s" 1 . To more clearly illustrate this point, the max-| 
imum and minimum values of the H + velocity in each bin^arC. 
shown. Note that while the minimum values arc about the safn<^ 
for all values of the 0 + velocity, the maximum H + speed increase^ 
dramatically as the 0 + velocity increases. At altitudes above j 
2000 km it is expected that the plasma is collisionless. Thus the*j 
processes which arc responsible for the correlations seen in Figure - 
7 must cither be taking place lower in the ionosphere or are affect- 
ing the accelerating, ambipolar electric field. 



: ig . 7 . The average velocity of H * as a function ot the O * velocity for all 
lata from altitudes greater than 2000 km. The inverted triangles above the 
iverages are the maximum H ~ velocity values in each O * velocity bin The 
riangles below the averages are the minimum H * velocity values in each 
/elocity bin. 
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Velocity theory comparison. The average ion velocities dis- 
cussed and shown above exhibit large statistical variations. It is 
important to understand whether these variations are indicative of 
experimental uncertainty or actual geophysical fluctuations. 
Figure 8 is a scatterplot of all the H + velocities used in this study 
as a function of altitude. Two points should be noted regarding 
these results. First, the values tend toward zero as the altitude 
decreases and the width of the spread in velocity increases with 
altitude. Theoretical results suggest this type of trend. For 
example, the results of Rain et al. [1975] show the upward H + 
velocity increasing from a value near zero at —800 km up to a 
maximum value, in the altitude range between 1000 and 2000 km, 
which depends on the boundary condition imposed at 3000 km. 
The solid curve in Figure 8 is their result for an upper boundary 
velocity of 20 km s -1 . This curve encompasses the RIMS results 
and follows closely the altitudinal profile suggesting that the fluc- 
tuations are mainly geophysical in nature. 



1 000 2000 3000 

Altitude (km) 


4000 


Fig, 8. Velocity of H*" as a function of altitude for all data. 


Average Density' 

The average density as a function of altitude (for altitude 
intervals of 500 km) is shown in Figure 9. In this altitude range O* 
is always the dominant ion being of the order of IO 3 cm' 3 , while 
H + is — 100 cm -1 and He + is —10 cm -1 . According to polar wind 
theory the major ion, 0 + , will establish the polarization electric 
field and the minor ions, H + and He + , should have the same alti- 
tude dependence as the dominant ion. This is indeed the case in 
the observations. When normalized to the 0 + density, the H + and 
He + densities at the higher three altitudes agree with each other 
and with the 0 + to within 10%. The ions have different altitude 
variations at lower altitudes because of differences in their ion 
chemistry. 

As with the average velocities, the densities have been aver- 
aged for the cases of high and low F, 0 7 flux and summer and 
winter season. In the case of 0 + (Figure 10) the density increased 
with increased F| 0 7 flux at all altitudes. This increase (by approxi- 
mately a factor of 2) is consistent with an increase in the solar 
EUV flux which is assumed to accompany a long-term increase in 
^10 7 flux. Conversely, there was no significant change in 
either the H + density (Figure lDorthe He + density (Figure 12). In 
ihe case of H + , while the 0 + density increases from solar mini- 



Alfitude (km) 


Fig. 9. The average density of each ion species as a function of altitude for 
all data samples. 



Altitude (km) 

Fig. 10. The average density of O '* as a function of altitude for F 7 < 140 
(I0-- W m~ 2 Hz“ l ) and f, 07 > 190 (IO' 22 W m~ 2 Hz 1 ). 
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Fig. II. The average density of H* as a function of altitude for Fun < 140 
(10 ” W nT- Hz 1 ) and F, u 7 > 190 (I0 “ W m : Hz '). 
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mum to solar maximum conditions, the [H ] decreases and the 
results is little variation in H+ density. Likewise, Rain et al. 
[19786] found no significant solar cycle effect in the He + density. 

The seasonal variations are shown in Figures 13, 14, and 15. 
The 0 + observations revealed a decrease in overall density from 
summer to winter of about a factor of 3. The winter density at the 
highest altitudes observed dropped to <200 cm' 3 which had signi- 
ficant impact on the ability of RIMS to accurately measure the 
density. This is evidenced in Figure 14 which shows a large dif- 
ference in the observed H* density near 4000 km. In this altitude 
bin the summer density is a factor of 4 higher than the winter 
value. Conversely, at altitudes near 2000 km the densities were 
either nearly equal or higher in the winter. A positive electrical 
potential on the spacecraft resulting from the relatively low 
electron density at 4000 km in the winter has screened out a por- 
tion of the H * . Thus these observations represent a lower limit for 
the actual values. Note that at the lowest altitudes the winter H + 
density is higher than the summer value. This is consistent with a 
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Fig. 13. The average density of O" as a function of altitude for the summer 
season and the winter season. 
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Fig. 14. The average density of H + as a function of altitude for the summer 
season and the winter season. 
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Fig. 15. The average density of He * as a function of altitude for the sumn 


season and the winter season. 


winter enhancement in the neutral hydrogen density. Similarly ^ 
Hc + , at the lowest altitudes the winter density is significant! 
higher than the summer density, indicative of the winter neuu. 
helium bulge. But. as with the H + , the observed He + wint 
density at higher altitudes falls to the level of the summer observ 
tions due to the positive potential. 

Summary 

The classical polar wind picture is one of steady outflow 
light ions through a stationary background of 0 + Results of th 
study show that when the real outflows are averaged over per 
of months such steady pictures do emerge. On the other hand, t 
lack of correlation between various geophysical paramet 
which vary over different time scales ranging from hours' 
months, suggests that much of the variability in the polar cap c 
flows is a result of short time scale (< hours) phenomena. Ne 
theless. there do exist some interesting aspects of these long- 
averages. Seasonal variations in the light ion outflows follow 
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trend expected from known variations in the neutral atmosphere 
Namely, the flux of He * increased by a factor of 10 from summer 
to winter. As evidenced in the observations of the density and 
velocity this increase was primarily due to increased density in 
winter, a direct result of the winter [He] enhancement. The neutral 
hydrogen density is likewise known to maximize during winter. 
However, the H* flux had a maximum value in the spring. This is 
probably the result of the interplay between the seasonal variation 
of [ H ] and that of O* which exhibits a summer maximum. 
Combining these two cycles can yield a spring maximum for the 
product of [H] and n 0 +. 

Correlations with geophysical parameters revealed some vari- 
ations though most were small. Variations in the flux with 
changes in the solar F i0 ? cm flux were less than a factor of 3 for 
H*. However, changes in the velocity were evident at middle 
altitudes (i.e., 2000-3000 km) where the light ions were observed 
to undergo the largest net acceleration. The solar output, as 
characterized by the F i0 7 value, influenced the altitude at which 
this acceleration began, the altitude being —500 km higher at high 
F10 7 compared to low F| 0 ?• This effect is a result of the higher 
neutral temperatures in the exosphere due to the increased solar 
output. These higher temperatures increase [O] and decrease [H] 
|cf. Richards and Torr, 1985]. The increased [O] and higher 
photoionization rate result in more 0 + as observed (see Figure 10). 
As a result H *-0* friction, which restricts the H * upward speed, 
becomes important at higher altitudes. 

Conversely the maximum speed for the light ions is observed to 
be more sensitive to the season than the F| 0 .7 flux. The season, or 
more importantly, the solar zenith angle, is important in determin- 
ing the ionization rate and hence the photoelectron escape rate and 
the magnitude of the polarization electric field. Higher ion and 
electron temperatures during summer may also play a role 
depending on whether or not the H* flow is supersonic. 

The ion outflow velocities revealed the influence of solar input 
as discussed previously. There was also observed a high degree of 
correlation between the velocities of the different ions. As 
expected the light ions were well correlated in the magnitude of 
their upward speed. The ratio of the He* velocity to that of H + 
was 0.28 (with a correlation coefficient of 0.71). If both ions were 
energized in the ambipolar field under conditions of zero colli- 
sions and no gravity, this ratio should be 0.5. However, because 
of the mass dependence of the energy loss to gravity the net 
energy gain for He* is less and the ratio is lower. 

The upward speed of H * was also observed to be correlated with 
the direction of O* flow. While the average H* velocity did not 
show substantial variation with O* velocity, the range of H* 
velocities increased as the O* velocity changed from downward to 
upward. That is, higher upward H * speeds occurred at times when 
O* was also moving upward. These relationships are indicative of 
the interplay, by means of the polarization field and collisions, 
among the ions which participate in the polar wind flow. For 
example, downgoing O* may occur in regions where there exists a 
time-dependent decrease in the electron and/or ion temperature of 
sufficient magnitude to affect the polarization field. Thus while 
the partial pressure of H + in the ionosphere might be sufficient to 
cause outflow the reduced ambipolar field would restrict the mag- 
nitude of that outflow. Such a case has been modeled by Schunk 
and Sojka [1989] where they found 0*/H* counterstreaming in 
regions of decreasing temperature. 

The question of subsonic versus supersonic flow is beyond the 
scope of the current work, depending on a thermal analysis of the 
outflows. Some inferences can be made, however, by comparing 
Ihe observations of the H* velocity as a function of altitude to 


theoretical results (e.g.. Banks and Holier , 1968]. Solutions to 
the equations for the H + Mach number show profiles, correspond- 
ing to subsonic solutions, which reach a maximum at some alti- 
tude and either decrease or remain constant above that altitude. If 
the results in Figure 3 are taken as representative of such a situ- 
ation then the implication is that on average the maximum flow 
speed during winter was less than the sound speed. The conclu- 
sion may then be drawn that at times during the period of the 
observations the plasma pressure remained finite along these polar 
cap field lines. Indeed, several of the velocity profiles from 
individual orbital passes through the polar cap show a similar 
structure. Attempts to correlate this effect with geophysical and 
solar parameters have, to this point, been unsuccessful. To date, 
no theoretical work has addressed this issue through systematic 
investigation of outflow characteristics as a function of the 
pressure at the “far' end of the flux tube. 

There remains much to leant regarding the nature of the polar 
wind. The question of the polar wind Mach number is left 
unanswered by this study due to the lack of knowledge of the ion 
temperature. Techniques for the determination of ion temperature 
from all available RIMS polar ionosphere observations are being 
developed. However, the implication of a subsonic polar wind 
(e g., nonzero “backpressure" along polar cap field lines) should 
provide impetus for a detailed parametric study of the character of 
the polar wind velocity under such conditions using the currently 
available theoretical models. Future observations should shed 
more light on this topic. Likewise, observations made on shorter 
timescales will aid investigators in addressing the influences of 
rapidly varying geophysical parameters. Such measurements will 
be available from instruments such as TIDE on the ISTP/POLAR 
mission. 
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Update on polar ’.vind thermal characterization scheme 

The application of thermal analysis techniques :o the DE 1/RIMS polar wind dataset 
is an interesting challenge. While the RIMS design was a significant improvement to 
previous Retarding Potential Analyzer (RPA) instruments, its mess-resolving capabilities 
introduced additional complexities in the form of an cncrgy-dcpcndcnt solid angle aperture 
and a finite energy bandwidth beyond the RPA setting. In addition, in-flight anomalies 
produced significant difficulties in comparison to the analysis of conventional RPA 
datasets. Particularly relevant to the RIMS polar wind dataset is that the radial head RPA 
mechanism failed to operate properly early in the mission lifetime. 

Fortunately though, by adopting a numerical quadrature approach to the solution of 
.the response integral, using sums over the aperture area, solid angle, and energy limits, 
energy-dependent complications may be sufficiently simulated. And, because separate spin 
and RPA profiles are available for each of the three instrument heads, the radial head spin 
curve and the end head RPA curves are usable for temperature analysis even in the absence 
of the radial head RPA information. Since the effect of the radial RPA failure is as if the 
retarding voltage was constantly zero, the total counts measurement remains accurate. 
Numerical algorithms of this type tend to be computationally inreisive, hut have the 
advantage of being applicable to a variety of plasma distribution functions. 

In the case of the polar wind dataset, further considerations on the data analysis 
technique are necessary. Strictly speaking, spacecraft spin angle distributions observed in 
situ are not true pitch angle distributions. The measured spin angle is defined in the 
tefercace frame of the spacecraft, which is different from the plasma rc3t frame due to 
spacecraft motion (2 lu 8 kilometers per second, depending on spacecraft altitude). This 
difference in reference frame would be uf liule consequence provided the thermal speed of 
the plasma is large compared to the spacecraft "ram speed. For a spueeoi a.\ median speed 
of 5 km/s, a factor-oi-five ratio of thermal speed to spacecraft speed implies 3.5 eV for H~, 

13 eV for He + and 52 eV for OT Because polar wind energies are considerably below this 
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threshold, the RIMS spin angle data contain information aoo it both the angular distribution 
of the plasma in its own frame of reference (the conventional pitch angle distribution) and 
the drift motion of the plasma relative to the spacecraft. These two aspects of instrument 
response tend to dominate independently - one at times obscuring the other - depending on 
the plasma's thermal speed and its relation tc the spacecraft and plasma drift speed. 

Roughly speaking, it is for subsonic (thermal speed greater than drift speed) ion 
Mach numbers that the angular characteristics ui the distribution aic best icaulved and fci 
supersonic Mach numbers that drift motion is most easily resolvable. Note that in the case 
of H+, flow will be subsonic for thermal components as low as l eV when drift velocities 
are less than 15 km/s. The bottom hire, for KIMS polar wine data, is that the relative drift 
motion will dominate the spin angle response, leading to a "rammed” flux peak in tne 
direction of relative wind as seen in the spacecraft frame. The angular (ar.d hence, 
temperature) information will be apparent only in the detailed shaoe of the rammed pealc 
and may, at times, become indeterminate. The end-head RPA profiles then, become 
increasingly important to the accurate determination of our polar wind thermal properties. 

With all these cautions in mind, a numerical simulation model of the RIMS 
response to a flowing plasma distribution was written [Giles et ai . , ' 993]. The simulation 
technique improves on previous work by realistically modeling the interaction of a moving 
spacecraft with a drifting plasma and by allowing insertion of a variety of distribution 
functions (for example, bi-Maxwchian, generalized Lorentz (Kappa), and Spiizsr-Hdrm 
[1953]), the tluee -component drift velocity, and a variable spacecraft potential. For the 
polar wind study, the distribuliuii function is assumed to be a flowing bi-Max weiiiar. which 
Is accelerated through the potential difference between the spacecraft and the plasma. The 
potential crop i 3 assumed to occur in an infinitesimally thin sheath about the spacecraft and 
the detector, rendering the detail of sheath structure uniuipuitam. The implementation 
allows simultaneous fitting of the radial sector spin curve data and Lite cud sec .or RPA 
(energy) curves. 
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Figure 1 shows simulated RTMS spin-angle and RPA responses to a bi-Max wellian 
plasma distribution, (solid lines) and the RIMS polar wind data (crosses) for 1 982/106 
22:17 UT. Note that the instrument response, for a given detector look direction, is to ions 
with velocities anti-parallel to the detector look direction Plasma distribution parameters 


for this fit arc: 

density 224 an 3 

T|j 0.36 eV 

Tj. 0.28 eV 

field-aligned velocity 4.14 km/s 

anti-sunward convection velocity 2 km/s 


Note that these parameters correspond to a cold subsonic outwara flow (ion Mach 
number -0.49). Early predictions indicated that the H + polar wind thermal distributions 
should be anisotropic, with T|j>Tj. for supersonic outflows and T|j<Tj.for subsonic flows. 
The present example however, fits in with work by Wu and Taieb [1993] showing that an 
upward heat flux may exist even for subsonic flow at low altitude in tire high-latitude 
ionosphere. 

Processing of the RTMS polar wind dataset using :he RIMS simulation code 
continues as of April, 1995. Several studies nave boon completed which relieve difficulties 
associated with fitting multivariate functions. Frequently it is a matter of art to decide 
which variables in a model should be fit from the dataset, and which should be held 
constant at fixed values, for example values predicted by theory or derived in previous 
analysis. In the case of the polar wire data sc;, variables potentially free for fitting to the 
data are the spacecraft potential, plasma density, the two-componcnt plasma temperature, 
and the three -component plasma bulk flow velocity. In the present implementation, the ion 
field-aligned velocity component and flux measurements (and so, the density, value) are 
held fixed. Their values were obtained [Chandler etui , 1991] by integrating the observed 
flux over the spin cycle of the satellite for a zero RPA voltage. Alsu, the spacecraft 
potential value is heid fixed during the first stage of the fitting process, an initial value 
being derived from an empirical relationship between the plasma density and spacecraft 
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potential [Comfort et cil., 1988], Values fit from the data set then are the two-component 
temperature and the remaining two components of the bulk flow velocity. The 
perpendicular temperature is fit by chi-square minimization to the end-head PvPA profiles. 
The parallel temperature is similarly fit to the radial head spin angle profile. The orientation 
of the spacecraft for the polar wind data set is such that the expected anti- sunward 
cun vec lion velocity is almost completely described by the y component of the bulk flew 
velocity vector. This is fit first by minimization tu the taJiul head spin angle piofile and 
then refined through minimization to the end-head profiles. Last, the spacecraft potential 
value is fit in the same order, to the same three profiles. The fit in figure 1 was obtained by 
this method. 
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